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PREFACE 


The purpose of this report is to 
present the accomplishments of the 
National Space Technology Labora- 
tories Earth Resources Laboratory’s 
R&T Program for Fiscal Year >983. 

The report includes program 
activities sponsored by the NASA 
Office of Space Science and Applica- 
tions and Spac< Station planning 
projects associate, with the Space 
Station Task Force. All Laboratory 
activities directly support overall 
NASA goals and objectives. 
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ABOUT THE COVER 

The lower Mississippi'Louisiana-Gulf of Mexico area <$ represent- 
ed in Landsat 4 Thematic Mapper sensor data acquired on 
September 16, 1982. The data have been subjected to principal 
component processing in order to create the high-contrast black 
and white image. The insert is an expanded data segment cover- 
ing the NSTL. The canal system associated with the Shuttle 
Engine Testing Complex is clearly delineated The brighter ele- 
ments in the upper portion of the insert indicate the recent con- 
struction area of the U S. Army Mississipoi Ammunition Plant 


CONTENTS 


INTRODUCTION 

SENSORS AND DATA SYSTEMS 

• Thermal Infrared Muitispectral Scanner 

• Advanced Sensor Design and Operation 

• Software Development and Management 

AGRICULTURE AND RESOURCES INVENTORY 
SURVEY THROUGH AEROSPACE REMOTE 
SENSING (AgRISTARS) 

» Large Area Land Cover Estimation 

• Land Cover Change Detection and Monitoring 

• Inventory and Monitoring Applications 

• Sensor Implementation and Evaluation 

• Conservation Practices Inventory 

• Soil Erosion Modeling 


• TEST AND EVALUATION 

• Corridor Analysis 

• Urban Studies 

• SPACE STATION SUPPORT ACTIVITIES 

• APPENDIX 

• Sources of Additional Information 

• ERL Publications 


APPLIED RESEARCH AND DATA ANALYSIS 

• Monitoring Semi-Arid Rangeland Degradation 

• Discrimination of Small Heterogeneous Mine Features 

• Soil Delineation Research 

• Archaeological Investigation 

• NSF Long-Term Ecological Research 

• Geological Mapping 

• Scene-to-Scene Registration 

• Topographic Corrections 

• TM Preprocessing Technique 

• Geobotanical Application 

• Scene-lo-Map Registration 

• Microwave Sensor Analysis 

• Mid/Thetmal IR Remote Sensing Data Analysis 

• Vegetation Stress Assessment 


JOINT RESEARCH PROJECTS 

• Crop Mensuration and Mapping 

• Timber Resource Inventory and Monitoring 

* Wetlands Productive Capacity Modeling 

* Crop Condition. Assessment and Monitoring 


INTRODUCTION 


The role of the National Space 
Technology Laboratories Earth 
Resources Laboratory (NSTL/ERL) 
encompasses: 

• Research and analysis leading 
to systematic methodologies 
for converting remote sensing 
digital data into useful, earth- 
related information. 

• Design, development, and 
operation of advanced air- 
borne sensors. 

• Development and utilization of 
state-of-the-art data processing 
hardware and software 
systems. 

This FY83 Research and Tech- 
nology Annual Report documents 
NSTL/ERI/s continued eminence in 
multisensor utilization, software 
development, and digital data anal- 
ysis. The numerous ERL investiga- 
tions reported on contribute to the 
overall NASA capability to acquire 
and use remotely sensed data in 
specialized areas of research. The 
joint research projects provide con- 
tinued interface with the non-NASA 
user communities in testing the 
technology for defined problems/ 
solutions. Space Station support ac- 
tivities provide an indication of 
NSTL/ERL functional flexibility. 
The report is divided into six major 
sections: 

• Sensors and Data Systems 

• Agriculture and Resources In- 
ventory Survey through Aero- 
space Remote Sensing 
(AgRISTARS) 

• Applied Research and Data 
Analysis 


• Joint Research Projects 

• Test and Evaluation 

• Space Station Support 
Activities 

Within NSTL/ERL’s realm of 
program activity. Fiscal Year 1983 
can be characterized as a year of 
transition in a number of ways: 

• A major thrust over the past 
decade for the NASA earth 
resources remote sensing pro- 
gram has been the develop- 
ment of analysis techniques 
for applying Landsat Multi- 
spectral Scanner (MSS) derived 
information to resource plan- 
ning and management. This 
era has witnessed a growth in 
the remote sensing community 
and a significant enhancement 
in the application of the tech- 
nology to resource manage- 
ment issues in the future. 
However, the earth resources 
program is at the edge of a 
quantum leap in capability 
with the availability of Land- 
sat Thematic Mapper (TM) 
data. The improvements in 
sensor performance, inherent 
in 30-meter spatial resolution, 
seven spectral bands (both the 
number and location in the 
spectrum), and eight-bit digital 
resolution, portend a new era 
of technology advancement, 
potentially as dramatic as the 
operational use of MSS data. 
The ERL is utilizing its oppor- 
tunity to put the new TM 
capability to the test, as evi- 
denced in portions of this 
report. 


• In the past, the earth re- 
sources program has been 
somewhat pragmatically 
oriented toward domestic ap- 
plications. Starting in FY83 
there has been a positive shift 
in emphasis with the develop- 
ment of a new programmatic 
theme focused on Land- 
Related Global Habitability 
Science Issues. The objective 
of this Office of Space Science 
and Applications (OSSA) pro- 
gram is to investigate long- 
term physical, chemical, and 
biological trends and changes 
in the earth’s environment, in- 
cluding its atmosphere, land 
masses, and oceans. The ERL 
recognizes the defined science 
issues and has set about re- 
orienting its R&D to be in 
concert with the emerging 
program. 

• In the past year, selected 
capability of NSTL/ERL has 
been brought to bear on 
NASA’s next potential major 
mission— -the Space Station. 
Participating in working 
groups and technology studies 
of the Space Station Task 
Force during FY83, NSTL 
plans continued support of 
these activities and a mean- 
ingful role in Space Station 
development once it becomes 
an approved program. 

In summary, FY83 has, to a 
large degree, been a year of change. 
As the contents of this Annual 
Report indicate, the NSTL. is in 
concert with these transitions. 
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SENSORS AND DATA SYSTEMS 


NSTL/ERL’s sensors and data 
systems development program incor- 
porates evolving technological ad- 
vancements into the design and 
development of advanced airborne 
sensors and interactive data analysis 
systems in support of overall NASA 
requirements and all research ac- 
tivities conduced by the Earth 
Rr,o ^rccs Laboratory. 

'iSTL /ERL’s activities in the 
development of interactive informa- 
tion system capabilities combine 
rapidly evolving hardware technol- 
ogy with state-of-the-art software 
systems to support both research ac- 
tivities conducted at NSTL/ERL 
and operational applications of 
remote sensing technology. Program 
emphasis during FY83 focused on 
upgrading hardware and software 
system capabilities for accommodat- 
ing TM and microwave data and 
multichannel classification oroce- 
durcs which take advantage of array 
processor technology. In addition, 
very low cost system configurations 
based on microprocessor technology 
were developed and checked in sup- 
port of multilevel information 
requirements. 

Tht following subsections high- 
light some of the more significant 
activities under the Sensors and 
Data Systems element of NSTL/ 
ERL. 


THERMAL INFRARED 
MULTISPECTRAL SCANNER 

FY83 marked the first full year 
under which the Thermal Infrared 
Multispectral Scanner (TIMS) per- 
formed in an “operational” status, 
although the late FY82 flight pro- 


gram and subsequent laboratory 
tests had revealed a few significant 
sensor problems. Consequently, 
even though FY82 had officially 
concluded the R&D phase of TIMS 
activities, FY83 saw an equal num- 
ber of laboratory tests and analyses 
that identified the two most impor- 
tant sources of data degradation in 
the system. 

TIMS is a six-channel thermal 
multispectral scanner capable of 
measuring target radiation in 400 
nanometer intervals from 8.2 
through 9.4 micrometers, and in 
800 and 1 ,000 nanometer intervals 
from 9.4 through 12.2 micrometers. 
Under laboratory conditions, noise 
equivalent temperature differentials 
of 0.05°C to 0.30°C are achievable. 
The system is operated in a Lear 23 
aircraft at altitudes ranging from 
6,500 feet to 40,000 feet above the 
target, producing ground resolutions 
of from 5 to 30 meters. 

TIMS’ uniqueness lies in its 
multispectral nature. Each of the 
six detectors measures thermal 
radiation, often as temperature in 
degrees Centigrade, of the target. 
Were thermal radiation the only 
contributor to the detected energy, 
then all six bands would produce 
the same results. However, the 
emissivity of the target being 
overflown is also a contributor to 
the energy being detected. Because 
the emissivity of any given object is 
nor constant across the 8.2-12.2 
micrometer range, slight variations 
in signal levels, primarily attribut- 
able to target emissivity, are in 
the output data. These small varia- 
tions allow investigators to deter- 
mine, by remote means, the appar- 
ent spectral emissivity differences of 


selected targets. 

A secondary benefit of TIMS 
results from its very high thermal 
sensitivity. High sensitivity was 
designed into the system because the 
contributions of emissivity are so 
small, and with a thermal sensitivity 
of 0.3®C or worse, apparent spec- 
tral emissivity detection would be 
virtually impossible. Consequently, 
the TIMS consistently produces sen- 
sitivities on the order of 0.1 °C or 
less. Thus, this allows single-band 
usage for detecting very subtle ther- 
mal variations. Detection of boat 
wakes, thermal effluents, thermal 
shadows, and archaeological phe- 
nomena are examples of this secon- 
dary aspect of TIMS. 

Laboratory tests during FY8j, 
in progress since system delivery 
in 1981, isolated the sources of data 
degradation caused by microphonic 
noise. The two areas slated for im- 
provements during FY84 are the 
preamplifiers and the optical 
system. 

The original design for the pre- 
amplifiers utilized a single circuit 
board, with six preamplifiers, 
located some 18 inches from the 
detectors. This circuit board also 
housed the circuits for providing the 
proper bias currents to eac'. detec- 
tor. The long cable runs, dictated 
primarily by a lack of space in the 
aircraft enclosure for a closer 
mounting location, produced ex- 
cessive electronic noise susceptibili- 
ty. A new preamplifier design, 
tested in late FY83, demonstrated 
marked improvement to noise rejec- 
tion, and will be incorporated in 
FY84. Rather than using a single 
circuit board design, the new effort 
will utilize discrete preamplifiers 
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located on the detector /dewar hom- 
ing itself. Sensitivity improvements 
of 2X to 3X are anticipated. 

The TIMS optical system is very 
large, with many optical compo- 
nents. Incoming energy must pass a 
rotating mirror, a focusing parab- 
ola, a 90 s reflector, a collimating 
parabola, a diffraction grating, and 
a very fast imaging lens before 
reaching the detector array. The ar- 
ray, at 0.01? inch in width, is situ- 
ated on a focal plane only 0.015 
inch wide. Movement of any single 
component in the optical path that 
will produce a change in excess of 
0.001 inch at the focal plane will 
result in a decrease in signal level. 
Normal aircraft low-frequency 
vibration can cause movement of 
the array, because of the required 
detector cooling dewar design, far 
in excess of this critical amount. 
Numerous system modifications, in- 
cluding structural strengthening of 
the dewar, were undertaken in FY- 


83. These modifications, when cou- 
pled with rigid control of aircraft 
speed and aerodynamic loading, 
have produced excellent quality 
flight data. Further modifications 
will take place in early FY84 to 
produce a final configuration which 
will produce satisfactory results in- 
dependent of aircraft flight induced 
vibration. 

TIMS had acquired data over 38 
test sites as of August 1983. During 
a 15-flight-day period in July and 
August, 28 successful missions were 
accomplished. Total coverage is rep- 
resented in Figure 1, which shows 
the locations of study sites covered 
during the current fiscal year. 


ADVANCED SENSOR DESIGN 
AND OPERATION 

During FY82, the Earth Resourc- 
es Laboratory undertook four ad- 
vanced sensor design studies, utiliz- 


ing Center Director’s discretionary 
funds. The four design studies were 
as follows: 

• Design of a Variable Resolu- 
tion Multilinear Array (MLA) 
Pushbroom Scanner Prototype 

• Design of an Active (Target- 
Illuminating) Multispectral 
Scanner 

• Design of a Scanning Terrain 
Profiler 

• Design of a Microprocessor- 
Controlled Field Spectrometer 

The Advanced Sensor Design 
activity for FY82 was formulated to 
take each study to a point which 
would allow the selection of one ac- 
tivity to be pursued, during FY83, 
to the point of production of proto- 
type hardware. 

The first design study was select- 
ed for the FY83 effort. As pror 
jected in late FY82, the purpose of 
this activity was to complete the 
design of a prototype MLA simula- 
tor and to fabricate and flight- 



Figure 1. FY83 TIMS Missions 
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evaluate the resulting sensor. The 
sensor was to be used to simulate 
the data anticipated from future 
satellite-borne MLA sensors, and to 
explore linear array technology as it 
applies to airborne remote sensing. 

Three spectral bands were select- 
ed to correspond to those projected 
for the future SPOT -satellite-borne 
instrument. Spectral bands, how- 
ever, were to be preflight selectable, 
and utilize thin-fllm bandpass Alters 
in screw-in mounts. A brord selec- 
tion of filters, with both very nar- 
row and very wide bandpasses, was 
envisioned. Overall detector spectral 
sensitivity would be limited to ap- 
proximately 500-1,050 nanometers. 

Early in FY83, it became obvious 
that the state-of-the-art for linear 
array charge-coupled devices was 
moving rapidly. The number of 
manufacturers of linear arrays had 
increased dramatically, with 1,728- 
element arrays being offered as pro- 
duction units. NSTL elected to pur- 
chase i 1 preproduction 3,456-element 
Texas Instruments array for testing. 
This unit demonstrated excellent 
characteristics under laboratory test 
conditions. With individual elements 
only 0.00042 inch square, the array 
appeared to be ideal for the vari- 
able resolution application. 

The prototype array was deliver- 
ed without thermoelectric cooling, 
so the initial testing was conducted 
only at room temperature. Signifi- 
cant variations in detector sensitivi- 
ty were noted as temperature 
changes took place, confirming not 
only the need for thermal control , 
but for accurate, thermally depen- 
dent calibration. 

Based on the complexity of the 
3,456-element array and the result- 
ing detail which must be applied to 
thermal testing, a concept of vari- 
able focal length foreoptics was 


dropped in favor of electronic con- 
trol over resolution. This meant 
that all data acquired from an air- 
craft altitude of 40,000 feet would 
produce ground resolutions of 5 
meters over an 8.7-mile-wide flight 
path. The sampling rate would be 
controlled on board, using additive 
buffers, to produce the required 
1Q-, 15-, 20- , 25-, and 30-meter 
ground cells. Operation at lower 
altitudes would still produce ground 
cells of approximately 1 meter. 
While this approach would com- 
plicate the on-board electronics, it 
would greatly simplify the optical 
design. 

By mid-August 1983, the design 
concept had produced the following 
definitions, cost projections, and 
schedule: 

• The prototype MLA device 
would consist of three 3,456- 
element linear array devices 
coupled to fixed focal length 
lenses to produce 5-meter data 
from 40,000 feet over an 8.7- 
mile swath, or 1 -meter data 
from 6,500 feet over a 1 .4- 
mile swath. 

• Spectral bands would coincide 
with those of SPOT, with al- 
ternate bands between 500 and 
1 ,050 nanometers available at 
the investigator’s discretion. 

• Projected hardware costs were 
S28.5K. 

In September 1983, NSTL was 
requested to investigate * change in 
scope of the overall MLA develop- 
ment project to determine its appli- 
cability to bidirectional reflectance 
measurements. Although modest 
off-nadir (fore/aft) viewing by one 
array for stereo effect had been 
designed into the original prototype, 
the bidirectional concept called for 
three channel simultaneous viewing 
of nadir and off-nadir angles rang- 


ing to 45* fore or aft. Additional 
complexities lie in the control of 
geometric distortions introduced by: 
(1) the large off-nadir viewing 
angle, and (2) the impact of aircraft 
motion on the data. 

A short, intensive study into the 
feasibility of fabricating a multi- 
linear array based bidirectional 
scanner and the associated dita 
management and utilization com- 
plexities was concluded by mid-Sep- 
tember, and produced the following 
definitions, cost projections, and 
schedules: 

• The Bidirectional Reflectance 
Imaging Array Radiometer 
(BRIAR) would consist of six 
3,456-element linear array 
devices coupled to fixed focal 
length lenses to produce 5- 
meter data from 40,000 feet 
over an 8.7-mile swath, or 
1-meter data from 6,500 feet 
over a 1.4-mile swath. 

• The six linear array lens units 
would consist of two spectrally 
matched pairs of three each. 
One of the pairs would be 
nadir looking, while the other 
set of three sensors would be 
capable of preflight selectable 
pointing angles of up to 45* 
either fore or aft of nadir. 

• Spectral bands would coincide 
with SPOT for both nadir and 
off-nadir pairs. Alternate 
bands between 500 and 1,050 
nanometers would be available 
at the investigator’s discretion. 
With both pairs pointing at 
nadir, six discrete spectral 
bands could be utilized. 

• Projected hardware costs were 


S53.0K, as follows: 


• Arrays 

$ 9. OK 

• Lenses 

12.0K 

• Circuit Boards 

12.0K 

• Mechanical Parts 

7.5K 
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• Thermoelectric 

Coolers 9.0K 

• Miscellaneous 3.5K 

• Flight tests were projected as 

being possible during the tenth 
month after adequate resourc- 
es are made available and the 
project is formally initiated. 
Should the NSTL Earth Resourc- 
es Laboratory elect to pursue the 
BRIAR approach, additional discre- 
tionary funds will be sought in 
FY84 to augment existing funds 
presently assigned to the prototype 
MLA effort. It is anticipated that 
additional Research and Technology 
Operating Plan (RTOP) funds 
would also be required to complete 
the bidirectional reflectance mea- 
surement task. Required manpower 
is currently being sought to comple- 
ment the existing design team. 


SOFTWARE DEVELOPMENT 
AND MANAGEMENT 

The Earth Resources Laboratory, 
over a number of years, has devel- 
oped a transferable set of software 
for processing imagery data from 
the Landsat Multispectra] Scanner, 
Seasat Synthetic Aperture Radar, 
airborne Thematic Mapper Simula- 
tor, and other remote sensors. 

These data can be combined with 
topographic, soils, rainfall, and 
other data to produce resource 
management and socioeconomic in- 
formation. ELAS, an acronym for 
Earth Resources Laboratory Ap- 
plications Software, is an operating 
subsystem that functions under a 
computer’s standard operating sys- 
tem. ELAS includes state-of-the-art 
capabilities for disciplinary user 
analysis, and was designed so that it 
can be transferred easily to most 
computers with 16- and 32-bit ar- 

• 


chitecture. Thus, ELAS uses a com- 
mon data file format with common 
system interface routines, making 
the addition of new program mod- 
ules relatively simple. 

ELAS is structured to provide a 
single point of contact between the 
computer and the user. As such, 
ELAS incorporates the necessary 
operating subsystem software with 
application software to accomplish 
a variety of tasks 'sociated with 
the processing of remotely sensed 
and other data types. The various 
modules address applications in 
agriculture, forestry, and range 
resources, geology studies, as well 
as land use. change detection, and 
wetland inventory assessment and 
monitoring projects. 

ELAS is used by discipline-ori- 
ented scientists and resource man- 
agers and is fuliy documented. It 
has been installed on most popular 
minicomputers used by universities 
and State and Federal agencies, and 
is in use in the private sector. The 
fundamental ELAS source code for 
several different brands of comput- 
ers (e.g., Perkin Elmer, Varian, 
SEL, VAX) has been made avail- 
able to the public through COS- 
MIC, NASA’s software dissemina- 
tion facility. The Perkin Elmer ver- 
sion of ELAS provides approxi- 
mately 200 separate executable pro- 
grams designed to be user-friendly, 
and is easily expanded to address a 
variety of remote sensing processes 
and applications. The processed 
data that are output may be dis- 
played on a line printer, CRT dis- 
play, electrostatic printer plotter, or 
color film product. 

Considerable effort has gone into 
the development of techniques to 
process selected remote sensor data 
through mathematical models. Sev- 
eral models have been developed, 


such as multitcmporal registration, 
computation of statistical signa- 
tures, maximum likelihood classifi- 
ers, georeference mapping, principal 
components, stepwise regression 
analyses, and terrain analyses. In 
addition, several models have been 
developed at ERL to assist the in- 
vestigators in applying wetland 
habitat technology (shoreline length, 
shoreline density, and distance to 
water), change detection, and spa- 
tial analysis techniques. 

The software capabilities of 
ELAS are ever-changing in order to 
meet NASA’s research and develop- 
ment requirements, and the capabil- 
ities will continue to expand due to 
ERL’s commitment to research and 
development activities. During 
FY83, 22 new ELAS modules were 
under development by the ERL. 

The module names and current 
status are listed below: 

• Modules currently under 

development: 

• TIRO — Reformats da;a 
from NOAA 6 and 7 
satellites. 

• GRAN — Gaussian random 
number generator. 

• Modules currently undergoing 

user-oriented validation: 

• TEXT — Performs texture 
analysis of data. 

• TDIM — Generates filter 
weights and their response. 

• FIL3 — Non-array proces- 
sor version of FILT, 

65 x 65 elements. 

• TERT — Calculates effec- 
tive radiant temperature 
from TIMS “housekeep- 
ing” data with no at- 
mospheric correction. 

• SPED — Performs spectral 
editing on polygons con- 
tained in PGF subfile. 

• Modules recently placed in 
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production use: 

• TOPO — Reformats NCIC 
tapes to EL AS format. 

• BLAP — Array processor 
version of FLAP. 65 x 65. 

• FLT5 — Frequency occur- 
rence of count values for 
channel of sensor data 
(IDI) and for individual 
cover type (ID2). 

• COMP — Compresses data 
from widely scattered areas 
of interest into the smallest 
possible data file. 

• TACT — Calculates target 
temperatures from TIMS 
data with atmospheric cor- 
rections applied. 

• TIMREF — Reformats 
decommutated TIMS data 
to ELAS format. 

• TCCT — Reformats Land- 
sat 4 Thematic Mapper data 
tapes (6250 or 1600 BPI) to 
ELAS data file format. 

• VPLT — ELAS plotting 
module for Versa tec Plot- 
ter. 

• SGDZ and SIDD — Digitiz- 
ing modules for generating 
point identification data 
(georef control points and 
polygon identifying points). 

• LOWTRAN — Calculates 
atmospheric transmittance 
and radiance for a given at- 


mospheric path based on 
user inputs such as altitude, 
pressure, temperature, etc. 

• MISP — Grids partially 
gridded or random input 
data. 

• COLOS — Locates 
minimum cost path between 
two cells in a study area. 

• CSCR1BE — Contours and 
plots data from MISP. 

• CSVF - Converts ELAS 
data files to single value file 
(SVF) for input to COLOS. 
Also converts output file to 
ELAS format for further 
processing. 

In addition, three modules— 
COMD, LABL, and PDDD— 
received significant upgrades during 
FY83. COMD and LABL now have 
the added capabilities for multi- 
image addressing, multi-image dis- 
play, true-color definition, and 
image-to-image copying. PDDD was 
modified to have plot axes inserted 
on frontal views, shading done only 
on axes in frontal views, and the 
vertical Z axis removed. These mod- 
ules are also currently in the valida- 
tion state. 

The modules described above 
were developed under a production 
software Configuration Control 
Procedure implemented in May 
1983 to ensure that changes in pro- 


duction software are made in an 
orderly manner. The procedure en- 
sures that the programmer and 
scientific user validate the software 
for technical accuracy, complete- 
ness, and documentation adequacy. 
Additionally, it requires ail software 
to be approved by the Configura- 
tion Control Board (CCB), which 
consists of six members representing 
software/hardware, scientific user, 
and contractor personnel. 

Three software libraries are 
maintained at NSTL/ERL: (i) 
Development and Modification 
Library, (2) Production Processing 
Library, and (3) Master File Li- 
brary. All software development ef- 
forts are contained in the Develop- 
ment and Modification Library. 

The Production Processing Library 
contains the modules necessary for 
all data processing activities. The 
Master File Library is a restricted 
account containing the master 
source files and companion flies for 
all modules used in either of the 
other two libraries. 

With the implementation of the 
Configuration Control Procedure, 
an environment for the development 
of software in an orderly and timely 
manner has been created, with the 
assurance that any modifications or 
additions to ELAS represent im- 
provements in the utility of the 
overall software package. 
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AGRICULTURE AND RESOURCES INVENTORY SURVEY 
THROUGH AEROSPACE REMOTE SENSING (AgRISTARS) 


AgRISTARS, a research program 
designed to benefit agricultural con- 
cerns by using aerospace tech- 
nology, is a joint program of 
NASA and the United States De- 
partments of Agriculture (USDA), 
Commc.ee, Interior, and State 
(Agency for International Develop- 
ment). The focus of the program is 
to develop satellite remote sensing 
techniques for a variety of agricul- 
tural applications and mapping of 
general land use and vegetation. 

NSTL/ERl participation in 
AgRISTARS began through two of 
the seven major AgRISTARS proj- 
ects. During FY83, all NSTL/ERl. 
AgRISTARS research tasks were 
consolidated under one major proj- 
ect entitled the Land Resources 
Research Project. In cooperation 
with the USDA Statistical Reporting 
Service (SRS) and Soil Comcrvation 


Service (SC'S), this project is com- 
prised .>f six specific research tasks: 

• 1 arge Area I and Cover Esti- 
mation 

• I and Cover ( hange Detection 
and Monitoring 

• inventory and Monitoring 
Applications 

• Sensor Implementation and 
Evaluation 

• Conserve' ion Practices Inven- 
tory 

• Soil Erosion Modeling 

LARGE AREA LAND COVER 
ESTIMATION 

The USDA Statistical Reporting 
Service regularly provides agricul- 
tural estimates of crops and live- 
stock production. They employ an 
elaborate ground sam, linj scheme 


to estimate agricultural production. 
Various investigations have reported 
improved crop estimates by inte- 
grating remotely sensed spectral 
data into this process. Additionally, 
th re is an advantage to producing 
estimates ol total land cover as well 
as cropland area from this routine 
operation for other agency informa- 
tion needs both within SRS and in 
USDA (e.g., SCS). The goal of this 
research is to improve these statist! 
cal sampling techniques and to in- 
crease the precision of full-frame, 
statewide estimates for the total 
range of land covers, while main 
taming the ful definition of sub- 
state and county-wide areas. 

Previous rest arch (NSTL/ERl 
Report l*>h AgRISTARS Report 
IX Yl (MOHR) under this task in- 
dicated that the spectral diversity 
among the land cover types was 
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under-represented by the study sam- 
ple, which represented 0.21b of the 
subject land area. Subsequently, 
research was conducted in Robeson, 
North Carolina, which utilized seg- 
ment data corresponding to about 
10% of the land area. This research 
also analyzed the results of using 
single-date versus multidate Land- 
sat MSS, and the results of stratify- 
ing by soils. The results of all four 
classification approaches (unitem- 
poral, multiteriporal, stratified, and 
unstratified) were presented at the 
December 1982 AgRISTARS techni- 
cal symposium (Stoner 1983). The 
analysis of the results revealed the 
stratified, multitemporal procedure 
to have the highest overall accuracy. 
The addition of a June Landsat 
MSS data set to an August data set 
had about the same effect on over- 
all classification accuracy as the 
stratification of the August data set 
by soils, although the unstratified, 
multidate approach was clearly bet- 
ter for crop discrimination. 

Another activity under this re- 
search task included the use of 
Landsat MSS data acquired ever 
Kansas to (1) expand the SRS 
methodology of crop area estima- 
tion to include forest, urban, range- 
land, and water, and (2) determine 
the utility of such land cover infor- 
mation to State agencies. The latter 
activity included the production of 
various map products for review by 
State agencies. The results were pre- 
sented at the AgRISTARS technical 
symposium (May, et ai., 1983). 

The collective experience gained 
from the North Carolina and Kan- 
sas studies will be applied to a 
statewide land cover classification 
for Missouri during FY84. During 
this activity. 67 additional segments 
from nonagricultural strata, mainly 


forestland, will be added to the 
regular SRS sample comprising 430 
segments. Multitemporal analysis 
will be conducted rather than the 
unitemporal analysis conducted with 
Kansas dam. Methods to better de- 
scribe the spectral diversity of land 
cover will be investigated. The use 
of a priori probabilities for pro- 
viding additional land cover infor- 
mation will be investigated and 
presented in the form of a cover- 
type odds map. This will provide a 
reliability index for each pixel. 


LAND COVER CHANGE 
DETECTION AND MONITOMNO 

Change detection investigations 
have been conducted with USDA 
SRS und'r the Land Resources Re- 
search Project of the AgRISTARS 
Prognim to develop methods for 
updating resource surveys as 
changes in land "over and use oc- 
cur. These methods were developed 
and tested with Landsat MSS data 
in an area of east-central Louisiana 
(Burns 1983) where forest land is 
being converted - cropland and 
pasture. The metnods involve co- 
registering multispectrai data sets 
for each study area into a multi- 
date, georeferenced data base. The 
radiometric differences between 
equivalent ground areas are then 
characterized in terms of land cover 
and/or land use change through 
various analytical techniques. These 
results are brought into the same 
data base along with digitized land 
cover change maps derived from in- 
terpretation of aerial photography 
acquired contemporaneously with 
each set of Ltndsat MSS data. The 
results from each change detection 
technique can then be evaluated 


relative to this ground reference to 
demonstrate the performance of 
each under the conditions present. 

During FY83 the methods were 
extended into Kansas where range- 
land i- being converted to irrigated 
cropland, and into Arizona where 
cropland is being converted to ur- 
ban residential areas. Results in 
Kansas indicate that within the one- 
degree by one-degree study area 
more than 60,000 acres of new 
cropland were established upon 
previous rangeland over a six-year 
period. Final results within the 
Arizona study area are expected in 
early FY84. 

Additionally, these experiences 
have suggested that area sampling 
frames could be updated by inte- 
grating the basic SRS geographic 
references unit (i.e., frame unit 
boundaries) with the most recent 
MSS-derived land use classification. 
This procedure eliminates the time- 
consuming process of reconstructing 
the entire sampling frame and was 
accomplished for the Louisiana and 
Kansas study areas. This computer- 
assisted land use stratification was 
demonstrated in conjunction with 
the operational update of the 
Arizona area frame for parts of 
Maricopa County in FY83 and is 
being evaluated by SRS personnel 
as a replacement for current 
methods. 

Other activities under this investi- 
gation in FY83 included evaluating 
Advanced Very High Resolution 
Radiometer (AVHRR) data from 
the NOAA-7 polar orbiting satellite 
for large area reconnaissance and 
targeting of problem areas. Data 
corresponding with the Louisiana 
verification data base indicate a 
high potential for discriminating 
and mapping the major land covers 
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in that area as well as flood plain 
and urban delineation using night- 
time therm* > data. In FV84, tests 
will include merging these data with 
Landsat 4 MSS and TM data in a 
resampled data base updated for a 
nine-year interval to compare the 
capabilities of each sensor for large 
area change detection surveys. 


INVENTORY AND MONITORING 
APPLICATIONS 

Missouri Land Use Inventory and 
Monitoring 

The objective of the Missouri 
Land Use Mapping Project is to use 
remotely sensed and ancillary data, 
in cooperation with the USDA Soil 
Conservation Service (SCS), to 
identify and measure the change in 
agricultural land use in areas im- 
peded by urban and suburban 
growth. Changes in land use over a 
10-year time span, July 1972 to 
August 1982, will be measured. The 
areas selected for the study are St. 
Charles and St. Louis Counties, 
Missouri, which encompass approx- 
imately 1,000 square miles. Studies 
in other pans of the United States 
have revealed that urban and subur- 
ban encroachment into rural areas 
and the subsequent land conversion 
have had a deleterious effect on 
agricultural systems and rural land 
use. Changes in agricultural land 
use will include the following 
studies: 

• The loss of prime and impor- 
tant agricultural land to other 
land use 

• The conversion of land use 
from wooded and pasture 
lands to cropland 


• A measure of potential soil 
erosion from agricultural land 
conversion 

In order to conduct these studies, 
a 10-year land use classification of 
the study area has been initiated 
using Landsat TM (August 1982) 
and MSS (August 1972 and 1982) 
data, aerial photography, ~nd an- 
cillary data. The ancillary data in- 
clude digitized soil maps, digital ter- 
rain tapes, digital Primary Sample 
Unit (PSU) data from the SCS Na- 
tional Resource Inventory (NRI), 
and ground verification data. 

Landsat TM data (August 1982) 
and NRI-PSU digital data have 
been processed to produce a land 
cover classification. Aerial photog- 
raphy (July 1982) of St. Charles 
County was obtained from that 
county’s SCS representative, and it 
has been processed for use in the 
digital image classification pro- 
cedures. Plate 1 shows land cover 
details depicted in the Landsat TM 
data as compared with older sources 
of land cover information. Landsat 
MSS scenes, aerial photography for 
St. Louis County, and ancillary 
data are presently being obtained 
for the project. 

Data processing during FY84 will 
entail geographic registration of 
sensor and ancillary data sets. Both 
supervised and unsupervised statisti- 
cal class development and maximum 
likelihood classification will be per- 
formed on multichannel data sets. 
Ground truthing procedures will be 
conducted to verify designated land 
use classes. A soil erosion data base 
will be created for the study area 
using data values for each of the 
Universal Soil Loss Equation 
fa 


Inventory and Mapping of 
Rangeland Brush 

The USDA Soil Conservation 
Service has identified rangeland 
brush inventory as one of its fore- 
most information needs which 
might be met using remote sensing 
techniques. This research task was 
undertaken to determine the poten- 
tial utility of Thematic Mapper data 
for determining brush type and den- 
sity, and to evaluate different tech- 
niques for extracting rangeland 
brush information from data ac- 
quired over New Mexico. 

A 6- by 35-mile area in Lea and 
Roosevelt Counties in southeastern 
New Mexico was selected by SCS 
personnel for intensive study. This 
area has considerable variety in 
brush type and density, with snake- 
weed, shinnery oak, mesquite, yuc- 
ca, sand sage, and catclaw being the 
important species. When it became 
apparent that Landsat TM data 
would not be available for New 
Mexico, Thematic Mapper Simula- 
tor (TMS) data were acquired on 
May 1, 1983, as t. substitute. 

Ground truth data were collect- 
ed from 105 sites during the first 
week of August. Canopy cover was 
estimated for each brush species 
and for grass. Analysis of these 
data will determine which band 
combinations and data transforma- 
tions will yield the greatest separa- 
bility of brush types and density 
classes. The benefits of including 
additional information on range site 
and spatial heterogeneity will be in- 
vestigated. The effect on separabili- 
ty by reducing resolution will also 
be determined. 
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SENSOR IMPLEMENTATION 
AND EVALUATION 

SAR Procedure Development 

During FY83 continued effort 
was placed on data processing and 
analysis to assess the utility of Syn- 
thetic Aperture Radar (SAR) data 
for crop and land cover estimation 
and mapping. X-band SAR data 
had been acquired on March 10, 
1981, for a truck garden area in 
Dade County, Florida; on June 1$, 
1981, for a rice paddy area in 
Acadia Parish, Louisiana; and on 
June 29, 1981, for a forest area in 
Kershaw County, South Carolina. 

The First steps in tne data proc- 
essing task consisted of a visual in- 
spection of the optically correlated 
image film, the digitization of im- 
age film, and the radiometric cor- 
rection of the digit’zed data to 
reduce speckle noise and across- 
track striping. The next step was to 
integrate (register or overlay) Land- 
sat MSS data and aircraft-acquired 
X-band SAR data with three polari- 
zations (HH, HV, and VV) to form 


a multisensor data set for analysis. 
Subsequently, spectral signatures 
were developed and the data classi- 
fied through multichannel pattern 
recognition using existing ELAS 
software, information synthesized 
from field observations was used to 
evaluate the classification. 

The investigation of the forestry- 
oriented Kershaw County data sets 
was documented in NSTL/ERL 
Report 213 (AgRISTARS Report 
DC-Y2-04374), “Analysis of Data 
Acquired by Synthetic Aperture 
Radar and Landsat Multispectra! 
Scanner over Kershaw County 
South Carolina During Summer 
Season” (Wu 1983c). The accuracy 
evaluation of the three classifica- 
tions using SAR, MSS, and SAR/ 
MSS data sets for the Kershaw 
County study area is given in Table 
I . The results were presented as a 
poster display at the AgRISTARS 
technical symposium in December 
1982. The results were also pre- 
sented as a part of the paper titled 
“Analysis of Synthetic Aperture 
Radar Data Over a Variety of Land 
Cover” at the 1983 International 


Geoscience and Remote Sensing 
Symposium (1GARSS ‘83) on Sep- 
tember 2, 1983, in San Francisco, 
California (Wu 1983a). 

The investigation covering Dade 
County, Florida, and Acadia Par- 
ish, Louisiana, data is documented 
in NSTL/ERL Report 218, “Analy- 
sis of Data Acquired by Synthetic 
Aperture Radar Over Dade County, 
Florida, and Acadia Parish, Loui- 
siana” (Wu 1983b). 

Also, in FY83 Shuttle Imaging 
Radar (S1R-A) data over Chicka- 
sha, Oklahoma, have been acquired 
and converted into computer-com- 
patible tapes. During FY84, the 
SIR- A data and Landsat 4 TM data 
will be employed to delineate sur- 
face features related to soil conser- 
vation practices. 


TM Procedure Development 

During FY83, analysis of the 
correlation between Thematic Map- 
per Simulator spectral response and 
percent closure in the forest canopy 
was completed for the San Juan 
National Forest, Colorado, study 
site. The analysis showed that the 
correlation coefficient was highest 
for the relationship between TMS 
Band 5 and canopy closure. Detail- 
ed results were presented at the 
December 1982 AgRISTARS techni- 
cal symposium, and in an AgRI- 
STARS Technical Report (Butera 
1983). 

Two TMS data sets of Kershaw 
County, South Carolina, were also 
analyzed at NSTL/ERL during 
FY83. The TMS data were collected 
during spring (May 1981) and win- 
ter (January 1982) seasons. The ob- 
jective of the project was to con- 
duct a supervised classification of 
each data set and to assess the ac- 


TABLE 1 

Kershaw County Land Cover Classification Accuracy 


LAND COVER CLASSES 


DATA SET USED 


MSS 4-BAND 3 POL X-BAND 
SAR OATA 

SAR a MSS 
BANDS 5, 7 

DECIDUOUS FOREST 

79 

79 

89 

PINE FOREST 

60 

44 

67 

PINE FOREST WITH SLASH 

65 

22 

68 

CLEARCUT 

63 

39 

79 

PASTURE/FALLOW FIELD 

64 

35 

77 

CROPLAND 

84 

59 

91 

WET FIELD 

70 

67 

86 

WATER 

96 

92 

99 

OVERALL ACCURACY 

72 

59 

82 
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Sf. Charles County, Missouri. Land use change along the Missouri River— 
woodland to cropland. 
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CIR Aerial Photography, 
September 1974 


Land sat 4 
Thematic Mapper Image 
August 1982 


HIGHWAY 40 


Sf. Louis County, Missouri. Land U3e change— woodland and cropland to 
suburban development. The completed Chesterfield Mall and new housing 
I development are seen along RoutJ 40 on the TM image. 
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curacy of TMS land cover classifi- 
cation with special attention to the 
accuracy of classifying various for- 
est types in both seasons. 

Field work was conducted in the 
spring of 1982 to describe homo- 
geneous land cover areas. These 
areas were located in the TMS data 
in order to generate training statis- 
tics which guided the supervised 
classification of both data sets. A 
separate set of verification poly- 
gons, reviewed in the field, will be 
used to assess the forest and land 
cover classification accuracy in early 
FY84. 

Thematic Mapper data collected 
during a November 16, 1982, Land- 
sat 4 pass were received in late 
FY83 at NSTL/ERL. In FY84, the 
TM data quality will first be evalu- 
ated to determine if the data (repre- 
senting the fall season) can be proc- 
essed and subsequently compared to 
land cover classification accuracies 
of the spring and winter TMS data. 

Submission of a final report ad- 
dressing the results of the seasonal 
land cover and forest classification 
accuracies is anticipated in Decem- 
ber 1983. 


CONSERVATION PRACTICES 
INVENTORY 

The goal of this research task is 
to determine the utility of remotely 
sensed data for the identification 
and inventory of existing soil con- 
servation practices. Five watersheds 
selected for the Rural Clean Water 
Program have been chosen as study 
sites in Mississippi, Oklahoma, 
Kansas, Illinois, and Idaho. Aerial 
photography (intermediate-scale and 
small-scale; and Landsat MSS data 
have been acquired for all five 
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study sites. One set of Landsat TM 
data has been acquired for each of 
the first four sites. SIR-A data have 
also been acquired for the Okla- 
homa study site. 

Preliminary data analysis during 
FY83 indicates that some conserva- 
tion practices can be successfully 
detected in Landsat 4 Thematic 
Mapper data (Plate 2). However, a 
number of existing practices are of 
such size and definition that present 
sensors and standard techniques 
cannot detect them with great ac- 
curacy. During FY84 special image 
enhancement techniques will be em- 
ployed with the hope of improving 
detection capabilities by highlighting 
patterns and properties of conserva- 


tion practices. Digitized soils and 
topographic data have been inte- 
grated with remotely sensed data 
bases, and will help focus attention 
on areas of critical interest. 

In addition, high-resolution pho- 
tography at the two different scales 
has been evaluated for the Missis- 
sippi, Oklahoma, and Kansas test 
sites. A matrix (Table 2) has been 
developed from the results of that 
interpretation which indicates the 
requirements for identifying a num- 
ber of different practices. Aerial 
photography from the Illinois and 
Idaho test sites will be evaluated for 
additional conservation practices in- 
dicative of those regions and results 
will be reported in FY84. 


TABLE 2 

Photo-Interpretive Matrix for Conservation Practices 
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SOIL EROSION MODELING 

The goal of this task is to survey 
the feasibility of using remote sens- 
ing techniques for providing various 
inputs to an erosion model using 
the Universal Soil Loss Equation 
(USLE). Four areas experiencing 
large soil loss rates have been cho- 
sen as study sites in Alabama, 
Mississippi, Oklahoma, and Kansas. 
Remotely sensed data from Landsat 
MSS and TM are being integrated 


with digitized soils and topographic 
data bases as input for determining 
the values of the various factors 
that make up the USLE. 

All remotely sensed data have 
been acquired and construction of 
digitized soils and topographic data 
bases has been completed in FY83. 
Preliminary use of the USLE model 
with Landsat MSS data in Missis- 
sippi (Plate 3) and Alabama have 
been reviewed with the USDA Soil 
Conservation Service personnel of 


those States and results indicate 
great promise. During FY84, results 
obtained from the model by using 
Landsat MSS data will be compared 
with results obtained when using 
TM data (or TMS as needed). A 
sensitivity analysis will be conducted 
to determine the influence and ac- 
curacy of the factors employed in 
the USLE model, whether they be 
derived from remotely sensed data, 
a soils data base, or a topographic 
data base. Results will be reported 
in FY84. 
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Conservation Practices Inventory 



Landsat Thematic Mapper scene 40136-16362 acquired November 29, 1982. Image is a 
computer-generated color composite using Bands 3, 4, and 5 to highlight soil conservation 
practices and erosional effects. 
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ORIGINAL PAS# 

COLOR PHOTOGRAPH 

AgRISTARS PLATE 3 

Soil Erosion Modeling 



IXIMUH SLOPE 


These images portray soil loss in tons/acre/year as determined by the Universal Soil Loss 
Equation (USLE) for the Goodwin Creek Watershed in Mississippi. The images use the max- 
imum and minimum slope values within soil mapping units to demonstrate the effect of slope 
on the USLE. In general, values greater than five tors/acre/year exceed tolerance limits. 


GOODWIN CREEK WATERSHED EROSION HAZARD 
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• APPLIED RESEARCH AND DATA ANALYSIS 

r 

• Monitoring Semi-Arid Rangeland Degradation 

• Discrimination ot Small Heterogeneous Mine 
Features 

• Soil Delineation Research 

• Archaeological Investigation 

• NSF Long-Term Ecological Research 

• Geological Mapping 

• Scene-To- Scene Registration 

• Topographic Corrections 

• TM Preprocessing Technique 

• Geobotamcal Application 

• Scene-To Map Registration 

• Microwave Sensor Analysis 

• Mid/Thermal IR Remote Sensing Data Analysis 

• Vegetation Stress Assessment 
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APPLIED RESEARCH AND DATA ANALYSIS 


NSTL/ERL’s Applied Research 
and Data Analysis activities concen- 
trate on research to understand the 
factors that affect or influence mea- 
surements made with remote sen- 
sors; i.e., the manner in which these 
factors manifest themselves in the 
data. During FY83, attention was 
focused on the middle and thermal 
infrared, and microwave regions of 
the spectrum. The Landsat Themat- 
ic Mapper (TM) and Shuttle Imag- 
ing Radar (SIR-A) sensors were 
drivers of this research, although 
some emphasis was placed on 
ground and airborne sensors such as 
the Barnes 12-1000 multiband field 
radiometer and the airborne The- 
matic Mapper Simulator (TMS) and 
Thermal Infrared Multispectra! 
Scanner (TIMS). 

Investigations were conducted 
with both fundamental and applied 
research orientations. Fundamental 
research tasks included studies of 
plant and soil temperatures, leaf 
water content and altering factors, 
vegetation stress, topographic cor- 
rections, atmospheric modeling, and 
data registration. Applied research 
tasks were directed toward soils 
mapping, arid vegetation monitor- 
ing, geological/geobotanical applica- 
tions, and archaeological investiga- 
tions. 


MONITORING SEMI-ARID 
RANGELAND DEGRADATION 

The goal of this research is to 
develop techniques for applying 
remote sensing to the assessment, 
monitoring, and prediction of semi- 


arid rangeland degradation. 

During FY83, field studies fund- 
ed by ERL were conducted by the 
Arid Lands Institute (University of 
Arizona) to determine degradation 
trends and to select vegetation pa- 
rameters that were both useful as 
indicators and potentially monitor- 
able by remote sensing. Results 
were reported in a paper entitled 
“Indicators of Rangeland Change 
and Their Potential for Remote 
Sensing” (P. L. Warren and C. F. 
Hutchinson, 1983 Journal of Arid 
Environments). 

The utility of MSS spectral in- 
dexes for monitoring rangeland 
vegetation was determined by ana- 
lyzing correlations between spectral 
indexes and critical vegetation pa- 
rameters and by determining rela- 
tionships between temporal changes 
in spectral indexes and changes in 
vegetation. MSS Band 3, albedo, 
and the Kauth-Thomas brightness 
component appear to be useful for 
monitoring total vegetation cover. 
Multiseasonal green vegetation in- 
dexes could be used to estimate 
change in the shrub/grass ratio. In 
the analysis of change in spectral 
indexes and vegetation over a four- 
year interval, spectral index change 
appeared to be offset from true 
change. Apparently the methods 
commonly used to standardize data 
sets for differences in solar eleva- 
tion effects and sensor radiometric 
response were not completely suc- 
cessful. Details of this study were 
presented in a report entitled 
‘‘Assessment of MSS Spectra! In- 
dexes for Monitoring Arid 
Rangeland” (Musick 1983). 


DISCRIMINATION OF SMALL 
HETEROGENEOUS MINE 
FEATURES 

Research continued during FY83 
on the extension of the thematic 
masking technique, developed in 
FY82 for the extraction of surface 
mine features in eastern Kentucky, 
from MSS data to TMS data col- 
lected for the same study area. A 
ratio of the first and second prin- 
cipal components derived from the 
data sets has been used to define 
spectral signatures exclusively for 
mined areas, with all other non- 
mine signatures ‘‘masked out.” 
These spectral signatures were then 
used to produce a thematic classifi- 
cation of surface mines within the 
eastern Kentucky study area. 

For an additional test of its 
spatial and temporal extension 
characteristics, the thematic mask- 
ing technique was applied to two 
dates of MSS data over a study 
area in the western Kentucky coal 
field. A thematic extraction of sur- 
face mines has been produced from 
April 1979 and May 1973 MSS 
data. The results show that surface 
mines within the western Kentucky 
study area have been successfully 
delineated and classified using the 
thematic masking procedure. 

This technique applied to extrac- 
tion of surface mine features was 
presented at the Association of 
American Geographers Annual 
Meeting held in Denver, Colorado, 
in April 1983. These final results 
will be formally documented in 
early FY84, which will complete this 
research investigation. 
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SOIL DELINEATION RE6EARCH 

The primary goal of this research 
is to develop remote sensing tech- 
niques capable of delineating soils 
in a manner which would serve to 
expedite the preparation of higher 
order soil surveys. 

Controlled ground-based experi- 
ments relating soil properties and 
field conditions to variations in soil 
spectral response form the basis for 
studies aimed at correlating remote- 
ly sensed data with soil map units 
and soil properties. Laboratory and 
in situ field measurements of reflec- 
tance are taken with a Barnes 
Model 12-1000 modular multiband 
radiometer, covering the six reflec- 
tive bands of the Landsat Thematic 
Mapper in addition to a near-IR 
band from 1.15-1.30 micrometers. 

Test sites were selected in arid, 
semi-arid, subhumid, and humid 
climatic regions to ensure that a 
variety of soil and climatic condi- 
tions would be represented (see 


Figure 2). Results from FY83 re- 
search at these sites are discussed 
below. 

New Mexico Site 

Sandy soils in this area are gen- 
erally more reddish than most fine- 
textured soils and are therefore sep- 
arable from most other soils by 
their spectral signatures in the TMS 
visible bands (Bands 1 , 2, and 3). 
However, reddish clay soils are <'cn- 
fused with reddish sands if only the 
visible bands are used (Section A of 
Plate 4). 

TMS Band 5 (1.55-1.75 microm- 
eters^ appears to be useful in sep- 
arating sandy from fine-textured 
soils. Compared to sands, flnc-tex- 
tured soils have generally lower 
reflectance in Band 5 relative to 
visible and near-IR bands. A com- 
bination of Band 5 with two visible 
bands thus separates soils by both 
texture and color (Section B of 
Plate 4). 



t - PAWNEE NATIONAL GRASSLANDS, COLORADO 

2 - JORNADA EXPERIMENTAL RANGE, NEW MEXICO 

3 - LITTLE WASHITA WATERSHED. OKLAHOMA 

4 - GENEVA COUNTY, ALABAMA 

Figure 2. Soil Delineation Research Study Sites 


in FY84 tests will be conducted 
to determine the utility of this band 
combination and others for soil 
mapping. Additional work will ex- 
amine the relationships between 
ground-based reflectance measure- 
ments and soil properties as deter- 
mined by laboratory analysis. Also, 
correlations between Thermal In- 
frared Multispectra! Scanner signa- 
tures and soils will be investigated. 

Colorado Site 

TMS, TIMS, and ground truth 
data were acquired for this site in 
late FY83. In this area of shortgrass 
steppe there is less exposure of bare 
soil than in the more arid New 
Mexico site. A major goal of in- 
vestigations at this site in FY84 will 
be to determine the extent to which 
partial vegetation cover alters the 
intrinsic spectral response of soils. 


Alabama Site 

Laboratory spectral response 
measurements of the surface and 
subsoil samples from four bench- 
mark soil series of the humid region 
indicated a high correlation between 
Band 7 (2.08-2.35 micrometers) and 
clay content (Figure 3). In situ field 
'Measurement of the same four soils 
series at the Alabama site were ob- 
tained at approximately the same 
time that TMS data were acquired 
over the area in Spring 1983 during 
a period of maximum bare soil ex- 
posure. A comparison will be made 
in FY84 of spectral responses from 
laboratory, field, and remotely 
sensed measurements. 

Oklahoma Site 

Digitization of the soils data base 
for Oklahoma was completed and 
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In these false color composites, Bands 1, 2, and 3 (Image A) and Bands 6, 2, and 3 (Image 
B) are simultaneously protected as blue, green, and red, respectively. Red clay and red sand 
have similar signatures in the 1/2/3 composite (A), but are readily separated in the 5/2/3 
composite (B) because Band 5 reflectance, relative to visible and near-IR bands, is generally 
high for sandy soils and low for clay soils. 

FALSE COLOR COMPOSITES FROM TMS DATA 
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Figure 3. Plot of Percent Ciey Content and Band 7 Spectral Response Values from 
Laboratory Data 


TM data acquired in FY83. Labora- 
tory and Held measurements from 
the Barnes radiometer are expected 
to be conducted in FY84, with em- 
phasis placed on the effect of crop 
residue in influencing soil spectral 
response. 


ARCHAEOLOGICAL 

INVESTIGATION 

In an era of dwindling budgets, 
remote sensing offers a rapid and 
inexpensive way to record and ana- 
lyze data previously unattainable 
through conventional archaeological 
survey methodologies. This type of 
data analysis is of major impor- 
tance in a scientific discipline which 
destroys its own data. Traditional 
methodology requires that explora- 
tory trenches be dug randomly 
within an area of interest, which 


basically results in a hit-or-miss 
process. Remote sensing improves 
this technique by rendering a 
t .seline understanding of the study 
area before expensive survey and 
excavation strategies are initiated. 
Although remote sensing has only 
recently been applied to archaeolo- 
gy, it promises to become a funda- 
mental part of the investigator’s 
scientific methodology in the 
1980’s, much as ceramic analysis 
and Carbon- 14 dating advanced ar- 
chaeological research in the 1950’s. 

During FY83, NSTL/ERL con- 
tinued its study to determine the 
feasibility of applying remotely 
sensed data from the Thematic 
Mapper Simulator and Thermal In- 
frared Multispectral Scanner to ar- 
chaeological applications. The 
Chaco Canyon region in north- 
west New Mexico was selected as 
the primary study area because it 


encompasses the largest concentra- 
tion of prehistoric ruins in North 
America. Particular emphasis was 
placed upon determining the loca- 
tion of prehistoric roadways, which 
are usually undetectable through 
standard archaeological survey tech- 
niques. These roads are generally 10 
meters wide and were constructed 
by the Anasazi culture around 900 
A.D. The reason for the existence 
of the roadways remains a mystery 
since the Anasazi did not possess 
the wheel nor beasts of burden. Ac- 
tive mineral exploration and devel- 
opr ent in this area of northwest 
New Mexico threaten to destroy the 
surviving remnants of the roadways 
within the next 5-10 years. 

Determining the location of the 
roadways is difficult from ground 
level. At dawn or dusk during cer- 
tain times of the year, some road- 
way segments stand out against the 
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landscape for a few brief moments. 
Despite the difficulty in viewing the 
roadways, their existence is never- 
theless supported through direct 
association with stairways, cause- 
ways, roadcuts, ramps, walls, curbs, 
changes in vegetation, ceramic arti- 
fact distribution, and alignment 
with outlying structures. Perhaps 
the major characteristic of the road- 
ways is their straightness of linearity 
regardless of the surrounding ter- 
rain. Through erosion, compaction, 
and vegetation differences, the 
roadways have, at best, survived as 
faint lines or segments of lines 
radiating out from Chaco Canyon 
across the surrounding desert. 

Image enhancement techniques 
using high-pass Gaussian filters 
were developed in order to accen- 
tuate the Chacoan roadway' on 
3-meter TIMS data. The Chacoan 
roads run predominantly in a north- 
south direction. To accentuate this 
phenomenon, directionally biased 
high -pass filters were used. These 
filters emphasized r.orth-south het- 
erogeneous features and de-empha- 
sized homogeneous features. A 
smoothing filter was subsequently 
used on the high-pass filtered data 
to remove noise and other non- 
associated roadway artifacts. The 
final product as presented in Plate 5 
reveals a well-denned Chacoan 
roadway network, although the 
roads are invisible to the human eye 
and are barely observable in the 
raw TIMS data. Results of the 
high-pass nitering techniques dra- 
matically demonstrated the ability 
of a directional filter to selectively 
remove noise along one axis of an 
image and extract desired informa- 
tion along the other axis. 

The development of filtering 
techniques offers unlimited potential 


in the analysis of prehistoric road- 
way and canal investigations around 
the world. High-pass filtering tech- 
niques applied to remotely sensed 
data can reveal the nature of trade 
networks and farming practices of 
ancient cultures in diverse environ- 
mental settings and consequently 
render a broader understanding of a 
civi‘:z;rion’s socioeconomic base. 
Recent discovery of Maya Canals 
and ridged fields in the Yucatan/ 
Peten have revised conventional in- 
terpretations of Maya subsistence 
patterns. Potential areas for appli- 
cation of linear analysis techniques 
include Inca roads in Peru, Roman 
roads in Palestine/Israel, ridged 
fields in Honduras, Maya roads in 
the Yucatan, and Megalithic roads 
in Britain. 

Documentation of the FY82 and 
FY83 research results will be avail- 
able in the first quarter FY84 as re- 
quired for NSTL Director’s Discre- 
tionary Fund reporting and will in- 
clude the work accomplished to 
date in Chaco Canyon, Poverty 
Point, Louisiana, and the Jackson 
Purchase area of Kentucky. Archae- 
ological investigations will continue 
during FY84 in cooperation with 
the American Schools of Oriental 
Research, Boston University, the 
American Institute of Archeology, 
National Park Service, and other 
prominent archaeology organiza- 
tions. 


NSF LONG-TERM ECOLOGICAL 
RESEARCH 

The National Science Foundation 
(NSF) funds research on long-term 
ecological phenomena at a national 
network of sites under the Long- 
Term F.cological Research (LTER) 


Program. LTER projects involve 
groups of investigators working at 
large, secure, and biologically di- 
verse sites. Systematic collection of 
comparative data on the biotic and 
abiotic components of these ecosy:- 
tems will support investigations of 
long-term hypotheses and questions 
in a series of core research topics. 

In FY83, two LTER projects 
were selected for coordinated re- 
mote sensing/ecological research: 

• Shortgrass Steppe: LTER. 

W. K. Laucnroth, et al., Col- 
orado State University. 

Site: Central Plains Experi- 
mental Range, Colorado. 

• interactions of Time and 

Space Variability in r. 
Chihuahuan Desert 
Ecosystem: LTER. W. G. 

Whitford, et al., New Mexico 
State University. Site: Jor- 
nada, New Mexico. 

In such projects, NSTL/ERL ac- 
quires remotely sensed data and 
develops techniques for extraction 
of information on vegetation and 
soil conditions. I TER investigators 
provide access to their data «md 
evaluate the extracted information 
with respect to their information 
needs. 

At the Central Plains site, color 
infrared photography acquired by 
NSTL/ERL is being used by LTER 
investigators in (heir intensive sur- 
vey of the surface geolorv, land 
forms, soils, and vegetation of the 
study area. Products from recently 
acquired TM and TIMS data will be 
used in this intensive survey, and 
variations in soil and vegetation 
parameters will be correlated with 
TMS and TIMS spectral signatures. 

TMS, TIMS, and color infrared 
aerial photography have been col- 
lected for the Jornada, New Mcx- 
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Pseudocolor composite image 
from TIMS • gh-pass filtered 
data, low-pass filtered data, and 
raw channel 3 (9. 0-9. 4 microm- 
eters) data, assigned to the 
blue, green, and red colors, 
respectively. The TIMS data 
were acquired at 5-meter reso- 
lution in August 1982 This im- 
age reveals prehistoric roads in 
Chaco Canyon, New Mexico, 
which were constructed by the 
Anasazi culture around 900 
A.D. Although these prehistoric 
roads are not visible from 
ground level, they demonstrate 
the ability of high-pass filters to 
extract information along one 
axis of an image using digital 
data. 




PREHISTORIC ROAD DETECTION 
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ico, site. The synoptic coverage of 
remotely sensed data will be used to 
extend the analysis of spatial and 
temporal variation from the inten- 
sive study site (a 3-km transect) to a 
larger area comprising the NMSU 
College Ranch, the Jornada Experi- 
mental Range, and adjacent private 
and public grazing lands. 

During FY84, in addition to con- 
tinuing research with NSF investi- 
gators in arid and semi-arid range- 
land ecosystems, cooperative re- 
search will be initiated in forest eco- 
systems with NSF-funded investiga- 
tors using the Andrews Experimen- 
tal Forest (Oregon), Coweeta Hy- 
drologic Laboratory (North Caro- 
lina), and Mt. St. Helens (Washing- 
ton) Long-Term Ecological Re- 
search study sites. 


GEOLOGICAL MAPPING 

This project, now in its second 
year, was designed to determine the 
utility of Thematic Mapper data (in- 
itially through simulation) in the de- 
tection and mapping of hydrother- 
mally altered rock that is character- 
istic of specific types of ore sys- 
tems. Emphasis is being placed on 
problems associated with extremely 
rugged terrain and varying degrees 
of vegetative cover. A second objec- 
tive is to develop the ability to in- 
tegrate the TM(S) data with other 
geologic data sets in a data base, 
interactively manipulate the data, 
and evaluate the results for their 
utility as a coherent whole. 

The study site chosen, in cooper- 
ate.. A'ith AMAX, Inc., as part of 
the NASA/AMAX research proje:t, 
was Mt. Emmons in Colorado. 

With more than 5,000 feet of relief 
in the area, and a very large, un- 


disturbed ore deposit, the study site 
is ideal for this project. 

In FY83, preliminary results of 
the project were published in Rick- 
man and Sadowski (1982) and 
Rickman (1983). These papers ad- 
dressed some of the practical prob- 
lems and coasideratin'.* in process- 
ing m vh da.a and some of the very 
preliminary products. Field work 
completed since the last yearly re- 
port strongly demonstrated the ex- 
tremely high overall quality and 
precision of the product images. 
Using these field observations as a 
base and guide, new images were 
generated to improve and supple- 
ment those previously created. 
Statistical analysis of the data has 
also revealed some apparent rela- 
tionship between vegetation density 
and several digital ratios. This rela- 
tionship was examined in detail in 
fourth quarter field work. If it can 
be oc‘ ,n tified, the relationship may 
be uscu to develop an ability to 
“see through” some vegetation to 
the rock and soil beneath. 

As this research effort nears 
completion, a review of results 
clearly shows the dramatic potential 
TM-type data have for those engag- 
ed in mineral exploration and min- 
ing. Using techniques detailed in 
Rickman (1983), a number of im- 
ages were generated from the same 
basic data, each fulfilling separate 
needs. For example, based on a 
priori knowledge of reflectance 
spectra, zones of anomalous hy- 
droxyl and ferric iron concentra- 
tion, free of vegetation, were 
detected by appropriate ratioing 
techniques, as shown in Section A 
of Plate 6. Such zones are known 
to correlate strongly with ore- 
bearing rock. A principal compo- 
nent composite (Section B of Plate 


6) emphasized zones that spatially 
correlate with known areas of 
mineralization. Other areas were 
also emphasized; field checking and 
a literature review have shown these 
to also be zones of mineralization. 

A vegetation map (Section C of 
Plate 6), developed from supervised 
statistics on ratios, functions both 
as a basic scientific tool and as a 
baseline environmental study. This 
image is more accurate and more 
detailed than the vegetation map fil- 
ed as part of the AMAX explora- 
tion Environmental Impact State- 
ment (EISj. Significantly, the image 
took a fraction of the time and ef- 
fort used to produce the map 
generated for the EIS. 

For the photo interpreter, the 
enhanced decorrelated composite 
(Section D of Plate 6) combines, in 
a near optimal presentation, infor- 
mation about all of the ground 
cover types. For the field worker, 
an excellent mapping base is provid- 
ed by the image in Section E of the 
plate. Here the geometrically cor- 
rected scanner data have been in- 
tegrated with digital elevation data 
expressed as contour lines, stream 
drainages, 7 '/: -minute quadrangle 
boundaries, and 2,500-meter tick 
marks based on the Universal 
Transverse Mercator (UTM) grid. 

Of general utility for a better 
understanding of an individual site, 
any of the products mentioned 
above can be presented in a per- 
spective view (Section F, Plate 6). 
This dramatically demonstrates the 
physical or topographic position of 
any area of interest. Such a presen- 
tation can significantly ease the 
problems of visualization. 

In FY84, with the expiration of 
the AMAX/NASA Letter of Agree- 
ment, the project will be concluded 
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with documentation of results and 
accomplishments in technical pa- 
pers. The TMS data set will be re- 
tained as a baseline or reference, 
but new processing is not antici- 
pated unless actual TM data be- 
come available. At that time, the 
process and techniques developed 
and used on the aircraft data will 
be applied to the satellite data, and 
the results will be compared. 

SCENE-TO-SCENE 

REGISTRATION 

In October 1981, NASA released 
an Applications Notice announcing 
the opportunity for participation in 
the Landsat D Image Data Quality 
Analysis Program. The objectives 
of the LIDQA Program were to 
quantify Landsat D sensor and 
systems performance from a user 
applications point of view. As its 
commitment to the LIDQA pro- 
gram, the Earth Resources Labora- 
tory has conducted scene-to-scene 


registration assessment research 
using Landsat 4 MSS data. 

The purpose of temporal registra- 
tion, using data from a sensor such 
as the Landsat MSS, is to produce 
a single multichannel data set from 
two or more individual multichan- 
nel data sets separated in time. The 
data from one date are used as a 
base and all other data are fitted or 
mapped onto them. It is generally 
accepted that by producing such a 
temporal data set, researchers can 
capitalize on changes which occur 
(or do not occur) over the time 
frame encompassed, thus gaining an 
added dimension that will lead to 
improvements in the utility of the 
data being examined. 

Results obtained from the use of 
temporal data are related to the 
degree to which two or more data 
sets can be registered one to an- 
other; the better the registration, 
the better the expected results, all 
other things being equal. Of par- 
ticular interest to this study were er- 
rors introduced by differences in the 


sensors being used; i.e., Landsats 1, 
2, and 3 versus 4. 

Findings from analysis conducted 
in this study are summarized in 
Table 3. The results obtained from 
the temporal registration of the 
1981 Landsat 2 MSS data to the 
base 1980 data set are typical of 
temporal registrations obtained in 
the past at ERL when using Land- 
sat 1, 2, or 3 data sets. The scan 
line and pixel (element) errors are 
quite similar, with no trends noticed 
in either magnitude or direction 
(based on a detailed analysis of 
overlay control points). Thus, a 
particularly uniform fit was 
achieved. 

The registration utilizing the 1982 
Landsat 4 MSS data (registered to 
the 1980 Landsat 2 data set) ex- 
hibited uncharacteristic error magni- 
tude in the pixel direction (when 
compared to both the corresponding 
scan line error of the same registra- 
tion and the pixel error of the 
Landsat-2-only registration). 

The source of this error was 


TABLE 3 

Results Obtained for the Registration of Three Scenes of Landsat 
P-Format MSS Data to a Common 1980 Landsat 2 Base Data Set 


COMPONENT 

1981 

LANDSAT 2 
(METERS) 

1982 

LANDSAT 4 
(METERS) 

1983 

LANDSAT 4 
(METERS) 

Pixel Error 

20.42 

40.32 

20.04 

Scan Line Error 

26.70 

26.45 

26.65 

Model RMS Error 

33.61 (74)* 

48.50 (65)* 

33.35 (72)* 


'Values in parentheses represent the number of points used to 
generate the corresponding statistic. 
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traced to scan mirror motion 
modeling coefficients in the God- 
dard Space Flight Center (GSFC) 
ground segment software. The 
problem was subsequently corrected 
at Goddard and the results are evi- 
dent in the 1983 Landsat 4 column 
of Table 3 (post-coefficient correc- 
tion), where the errors are in line 
with the Landsat-2-only values. 

It is concluded, therefore, that 
temporal registration of Landsat 4 
MSS P-format digital data with 
historical data obtained by Landsat 
2 can be performed to a degree of 
precision equal to that obtained 
when Landsat 2 data are used as 
both components of the temporal 
registration. 

Details of this study are con- 
tained in Anderson (1983). 


TOPOGRAPHIC CORRECTIONS 

Of necessity, the processing of 
remotely sensed imagery requires 
some assumptions about the data. 
Very important among these as- 
sumptions is that of uniform illumi- 
nation— the belief that all targets 
receive exactly the same irradiance. 
In any but the flattest of terrains, 
this is a very dubious assumption. 
Indeed, the amount of illumination 
received directly from the sun can 
vary by a factor of two in moder- 
ately hilly areas. Variation of this 
magnitude has obvious impact on 
the ability to successfully classify or 
analyze the data. For this reason 
ERL has engaged in a project to 
define the magnitude of the prob- 
lem, and adapt or develop practical 
techniques to correct the digital 
values in the imagery. 

An extensive literature review 
and analysis conducted in FY83 


revealed that the problem could be 
broken down into three major con- 
stituents: (1) the contribution to the 
recorded digital values due to atmo- 
spheric haze (the additive path radi- 
ance), (2) the variation of direct 
solar irradiance with slope and 
aspect, and (3) the percentage of 
the total irradiance from the sky 
hemisphere. 

From this understanding it was 
necessary to develop appropriate 
corrective software and evaluate its 
impact on the data. For this pur- 
pose, a data set covering the Lolo 
National Forest in western Montana 
was obtained. This package is an 
integrated data base including MSS 
data, 1 :24, 000-scale digital terrain 
information, and extremely detailed 
surface cover observation data. 
Using this data set and capabilities 
inherent in ELAS, the radiative 
transfer process was modeled and 
the raw MSS data (Plate 7, Section 
A) was modified based on the three 
constituents’ effects. This resulted 
in the image shown in Section B of 
the plate, which appears much 
“flatter” to the observer. The 


results of a more rigorous evalua- 
tion of the correction’s effect can 
be seen in Figure 4 and Table 4. In 
principle, after correction is applied 
to a given, specific land cover type, 
the digital values in an image 
should be independent of (lack cor- 
relation with) the angle of incident 
illumination (i.e., the cos i effect). 
Loio ground truth identified one 
representative class as being 95% 
tree and 5% shrub. Before correc- 
tion, the digital values for this class 
are strongly correlated to cos i 
(have a significant slope and R 
value). In other words, the mea- 
sured reflectance for a homogene- 
ous reflecting cover type appears to 
be varying with the angle of illumi- 
nation. After correction the cor- 
relation is almost totally gone, as 
evidenced by the R values (a mea- 
sure of correlation) in Table 4 and 
the horizontal shifting of the origi- 
nally sloped scattergram in Figure 
4. 

These and similar results suggest 
that the corrective processing is 
working appropriately. Further 
work is planned, both to further 


TABLE 4 

Lolo National Forest Ground Truth Sample 
Before and After Correction tor Topographic Effects 
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document the incremental value of 
each step in the formative correc- 
tion process and to study those 
results that appear anomalous. 

TM PREPROCESSING 
TECHNIQUE 

Extensive research has been con- 
ducted dealing with the application 
of feature selection techniques de- 
signed to facilitate the processing of 
Landsat MSS data. With the advent 
of the Thematic Mapper on Landsat 
4, the implementation of multitem- 
poral land use classification schemes 
or the use of band ratioing tech- 
niques to enhance the information 
in a unitemporal classification will 
stretch the capabilities of many ex- 
isting data processing systems. This 
will necessitate the development of 
feature selection techniques which 
reduce the dimensionality (number 
of channels) of TM data, while at 
the same time preserving the infor- 


mation content inherent in the raw 
data. Subset selection and feature 
extraction represent two different 
approaches to the dimensionality 
reduction problem. 

Subset selection employs multi- 
variate techniques to select the sub- 
set of original channels which pro- 
duce the best classification of a 
given data set. Feature extraction 
employs multivariate statistical 
methods to transform the raw data 
into new “axes” that are independ- 
ent of one another (orthogonality 
condition) and to classify the 
transformed data set based on the 
transformed values. The subspace 
of reduced dimension produced by 
either of these techniques is the 
feature space. By statistically com- 
paring the accuracy of the classifi- 
cation produced from the original 
measurement space with that gener- 
ated from the feature space, it is 
possible to choose the method that 
minimizes the probability of mis- 
classification. The other factor of 


concern is the cost (Central Process- 
ing Unit time utilized) inherent in 
classifying the feature space versus 
that incurred for the original mea- 
surement space. 

This study initiated in FY83 
employs multispectral scanner data 
from the Landsat 4 Thematic Map- 
per collected over a study site that 
surrounds the National Space 
Technology Laboratories in Han- 
cock County, Mississippi. The study 
site includes urban and residential 
areas, agriculture and pasture, for- 
ests, various types of fresh and 
brackish water, and a variety of 
wetland vegetation communities. 

TM data were collected in October 
1982. Extensive ground truth data 
have been gathered over this study 
site so the distribution of land use 
categories is well understood. 

The experimental design requires 
the comparison of the land use 
classification produced from six 
channels (visible and mid-IR) of 
raw TM data with the classification 
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Figure 4. Cosine of Incident Angle of Illumination (Band 6) 




After correction for 
atmospheric haze and 
angle of incident 
illumination 
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generated by processing the reduced 
number of channels resulting from 
the application of different tech- 
niques of feature selection. Subset 
selection methods, which allow the 
investigator to choose those chan- 
nels required to delineate the land 
use classes of interest, will be ex- 
amined, based on the raw TM data 
that have been grouped into train- 
ing statistics for known land use 
classes. Feature extraction methods 
will also be used to transform the 
raw TM data into new variables in 
a reduced feature space, which will 
then be classified utilizing the tech- 
niques available at the Earth Re- 
sources Laboratory. 

FY83 accomplishments include 
the following: 

• Receipt in March 1983 of a 
TM data set suitable for proc- 
essing (data collected in Oc- 
tober 1982). 

• Location and spectral editing 
of more than 160 ground truth 
sites in the TM data set, rep- 
resenting all the major land 
cover types of interest in the 
study. Spectral editing was 
necessary to ensure spectral 
homogeneity of land cover 
within the polygons used to 
delineate the ground truth 
sites. 

• Classification of the raw six- 
channel data (the thermal 
channel was not used), and in- 
itiation of accuracy verifica- 
tion of the resultant product. 

• Generation of a principal com- 
ponents transformed data set 
(six axes) and initiation of 
ground truth site integration. 

Anticipated FY84 activities include: 

• Completion of the evaluation 
of the six-channel raw data 
classification. 


• Completion of the principal 
components analysis. 

• Initiation and completion of 
canonical analysis and subset 
selection. 

• A comparison of all three 
approaches 


GEOBOTANICAL APPLICATION 

This investigation, under FY83 
RTOP 677-42-04, was initiated to 
develop and evaluate techniques for 
using the Thematic Mapper (initially 
through simulation) and other air- 
and space-borne systems for geobo- 
tanical mapping. The emphasis is 
on ore-bearing terrains in areas 
which are moderately to heavily 
vegetated. Specifically, emphasis is 
placed on differences in the vege- 
tative cover that appear to relate to 
the surficial geology. These varia- 
tions can be either differing vegeta- 
tion types due to geologic structure 
or perhaps subtle differences in 
canopy spectral reflectance caused 
by varying mineralization in the 
soil. 

The area of investigation centers 
around the Haile gold mine which 
is located near Kershaw, South 
Carolina. Thematic Mapper Simu- 
late data were acquired over this 
area in April and October 1982. 
Preliminary processing of these data 
removed noise and bad scan lines 
and corrected for overscan. The 
data were then georeferenced so 
that each raw data channel could be 
assigned to a channel in a geo- 
graphic information system. Geo- 
logic map data were also digitized 
and included as a channel of data. 
Processing the data using a classifi- 
er whose output simulates a three- 
band color composite, where the 


digitized intensity of each assigned 
color is dependent on the intensity 
values in each channel, yields a 
false color image with 236 distinct 
classes. An example of this is 
shown in Plate 8. 

The plate consists of the three- 
color composite image from the Oc- 
tober TMS data in which the data 
channel containing the digitized 
geologic map data has also been 
superimposed. In the digitized geo- 
logic data, divisions between major 
formations are shown in yellow, 
while narrow linear features, such 
as diabase dikes, are shown in 
white. 

In this plate, four major litho- 
logic divisions can be seen. The red- 
dish area in the lower right quad- 
rant corresponds closely to the Cre- 
taceous sediments (unconsolidated 
quartz sand), whose boundary with 
the underlying felsic metavolcanic 
rocks is shown in yellow. Note that 
the areas of felsic metavolcanic 
rocks are well depicted by the 
greenish areas of this image. The 
darker areas in the upper right 
quadrant correspond well to the 
prophyritic, coarse-grained adamel- 
lite of the Carboniferous Pageland 
pluton, while in the extreme upper 
left quadrant, the medium-grained, 
foliated, and metamorphosed me- 
tatonalite is shown in darker green 
and darker reddish coloration. This 
area is also separated from the 
felsic metavolcanic rocks by the 
yellow boundary. Diabase dikes are 
depicted in the geologic overlay as 
the linear northwest-southeast 
trending features shown in white. In 
the original data, a number of such 
features were visible (usually as 
some shade of green) but here are 
overlaid by the geologic graphic. 

FY84 studies will include the 
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analysis of products of data en- 
hancement techniques which have 
been applied to this area, utilizing 
specialized filtering techniques 
which will directionally enhance 
linear features. Additionally, similar 
data processing and analysis efforts 
will be concentrated on the northern 
portion of the study area (Uwharrie 
National Forest) where there is 
much less disturbance to the natural 
vegetation as a result of man's 
activities. 


SCENE-TO-MAP REGISTRATION 

The objective of this Fundamen- 
tal Research project is to assess the 
geometric accuracy of Landsat MSS 
and TM sensor products. A second 
facet of the project is to develop 
improved procedures for the regis- 
tration and rectification of Landsat 
data. 

An evaluation of the accuracy of 
U.S. Geological Survey (USGS) 

ER OS-generated P-format (geomet- 
rically corrected) data for Landsat 
MSS scenes from Kansas and Loui- 
siana was conducted in FY83 using 
the approach of Graham and 
Lucbbe (NSTL/ERL Report 197/ 
AgRlSTARS Report DC-Y1-04069, 
1981). The row bias was 251.8 pix- 
els for the Kansas data and -219.4 
for the Louisiana data, while the 
column bias was 100.3 pixels for 
Kansas and -95.6 for Louisiana. 
These results suggest that the P- 
format data for Kansas and Loui- 
siana are not accurately georegis- 
tered. Georegistration of P-format 
data utilizing Earth Resources Lab- 
oratory Applications Software 
(ELAS) for whole Landsat scenes 
utilizing 32 ground control points 
(GCPs) reported row bias values of 
0.03 pixel for Kansas and 0.17 for 


Louisiana; the column bias was 
-0.07 pixel for Kansas and -0.14 
for Louisiana. A similar study of 
A- format (radiometric corrections 
only) data, georegistered using 
ELAS software, gave a row bias for 
32 GCPs of 0.09 pixel for Kansas 
and -0.09 for Louisiana, with a col- 
umn bias of -0.15 pixel for Kansas 
and 0.44 for Louisiana. This work 
suggests that ELAS software pro- 
vides a more accurate georegistra- 
tion product than does the EROS 
P-format-generated product. This 
work is reported in the Proceedings 
of the NASA Symposium on 
Mathematical Pattern Recognition 
and Image Analysis (Dow 1983). 

Current work involves empirical 
studies to evaluate the influence of 
the number and spatial distribution 
(random, uniform, or clustered) of 
GCPs on the geometric accuracy of 
the scene-to-map registration proc- 
ess for selected portions of an 
A-format Landsat MSS scene. The 
A-format scene is georegistered 
using ground control points; the 
accuracy of the scene-to-map regis- 
tration process is evaluated by com 
paring the location of known 
ground reference points (GRPs) 
with the projected location of the 
points using the mapping equation 
produced in the georegistration 
process. The accuracy of this 
analysis is expressed in terms of the 
bias and standard deviation descrip- 
tors of Graham and Luebbe (1981). 
The goal of this work is to deter- 
mine the optimum number of points 
and spatial distribution of GCPs re- 
quired to georegister a given pro- 
portion of a Landsat scene. 

Future work includes analyzing 
the scene-tc ap accuracy of TM 
data from the R&D phase of the 
NASA LIDQA program. This 
analysis would proceed with an ex- 


perimental design similar to that 
employed with the P-format Land- 
sat MSS data. An artificial data set 
has been obtained from Purdue 
University which can be employed 
to evaluate the influence of GCPs’ 
number and spatial distribution on 
scene-to-map registration accuracy. 
This study design would be similar 
to that employed with the A-format 
Landsat MSS data analysis. 


MICROWAVE SENSOR 
ANALYSIS 

Under an Applied Research and 
Data Analysis (ARDA) RTOP task 
addressing Multisensor Technique 
Development, radar land cover 
analysis research objectives empha- 
sized the understanding of micro- 
wave attributes. The potential utility 
of radar imagery for surface feature 
identification and mapping estab- 
lishes a need for a basic under- 
standing of the radar signatures for 
various surface cover types and the 
use of this knowledge in the devel- 
opment of techniques oriented 
toward meaningful terrestrial infor- 
mation extraction. ERL obtained 
ascending and descending orbit pass 
data acquired over Baldwin County, 
Alabama, by the Shuttle Imaging 
Radar (S1R-A, L-band) and analyz- 
ed the data in conjunction with 
ground truth information to deter- 
mine which basic land cover physi- 
cal and/or biological properties af- 
fect microwave backscatter. S1R-A 
data were also integrated with 
Landsat MSS data to quantify the 
improvements in using a multisen- 
sor data set. 

Analysis of results indicates that 
cropland/pasture and water classes 
contain low digital count levels 
(related to radar backscatter), while 
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mixed pine/hardwood and urban/ 
inert classes contain high count 
levels for SIR-A data. The count 
levels of three pine forest classes are 
highly correlated with tree height 
and canopy structures. The use of 
MSS data improved separation of 
cropland/pasture classes from water 
and pine regeneration classes. The 
results of dassi fixation accuracy 
using SIR-A, MSS, and SIR-A/ 

MSS data are given in Table 5. 
These results arc included as part of 
a paper entitled “Analysis of SAR 
Data Over a Variety of Land 
Cover," published in the 1983 In- 
ternational Geoscience and Remote 
Sensing Symposium (IGARSS *83) 
Digest (Wu 1983a). 

After the launch of Landsat 4, 
data acquired by the Thematic 
Mapper over the Baldwin County 
study area in October 1982 were co- 
registered with SIR-A data to form 
a SIR-A/TM data set. The process- 
ing and analysis of these data to 
evaluate the complementary role of 
TM data in delineating forest and 
urban cover types will be continued 
and completed in FY84. In Sep- 
tember 1983, a data acquisition air- 
craft mission using the Jet Propul- 
sion Laboratory (JPL) L-band 
multipolarization SAR system was 
conducted. Analysis of these data 
from the agriculture- and forestry- 
oriented study area will be the ma- 
jor thrust in FY84. 


MID/THERMAL IR REMOTE 
8EN8INQ DATA ANALYSIS 

In order to effectively utilize in- 
formation derived from the mid- to 
thermal infrared regions of the elec- 
tromagnetic spectrum (0.76-14.0 
micrometers), research needs to be 
conducted dealing with the basic 


factors associated with land covers 
that influence mid/thermal IR re- 
flectance and emissivity. Fundamen- 
tal to such research is an under- 
standing of the relationships be- 
tween the life processes or physical 
state of a land cover and the phe- 
nomena that are measured in the 
mid/thermal IR region. 

Two areas of interest are being 
investigated in this ARDA RTOP 
project. The first deals with thermal 
studies (3.0-14.0 micrometers) and 
will focus on the estimation of 
plant temperature, which is not an 
end in itself but rather an initial 
point from which other studies re- 
quiring a knowledge of temperature 
can begin. The second area of in- 
terest, based at least in part on a 
target's physiological or chemical 
response to radiation with wave- 
lengths between 0.76 and 3.0 mi- 
crometers, is the determination of 
relationships that exist between leaf 
water content and data recorded by 
Earth-orbiting satellites. 


It has been demonstrated by 
numerous researchers that the 
potential exists to extract valuable 
information from data collected in 
the mid-IR region for plant/envi- 
ronment interactions, including re- 
search dealing with the reflectance 
changes due to leaf maturation and 
senescence, leaf type (gymosperm 
versus angiosperm), plant nutrition, 
site salinity, and water stress. In 
summary, the major influence on 
mid-IR reflectance by plant leaves is 
tied to the leaf water status, and 
any organism or condition which 
alters this status will have an effect 
on the reflectance as measured by a 
mid-IR-sensitive device. 

Interest in the thermal IR region 
of the electromagnetic spectrum is 
founded in part on the relationship 
that a target’s temperature has with 
its physical or physiological state, 
and how well the temperature can 
be used as an indicator of that 
state. The effect of temperatures 
can be related to the state of agro- 


TABLE 5 

Baldwin County Land Cover Classification Accuracy 


LANDSAT 
COVER CLASSES 

SIR-A 

DATA 

LANDSAT 
MSS DATA 

SIR-A AND 
MSS DATA 

Urban/Inert 

6.7 

35.7 

72.5 

Mixed Pine Hardwood 

77.9 

24.5 

83.3 

Native Pine 

91.0 

78.4 

83.0 

Pine Regeneration 

66.4 

41 .0 

78.8 

Cropland/Pasture 

74.4 

82.8 

93.1 

Water 

64.9 

98.5 

99.8 

Overall 

67.0 

68.0 

89.2 
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nomic crops: soil surface tempera- 
tures reaching 50-60* ’ 'an kill 
emerging plants; soil surface tem- 
peratures collected in the winter and 
early spring months may indicate 
the need for freeze protection prac- 
tices; and plant temperature in- 
creases with decreasing soil moisture 
availability and canopy tempera- 
tures can be used as an indicator of 
the relative need for irrigation. 
Temperature is also a useful quanti- 
ty to know in other fields, such as 
thermal discharge studies, soil- 
related investigations, geology, 
oceanography, geomorphology, at- 
mospherics, snow field studies, and 
urban studies. 

In the past, temperatures were 
obtained by making direct measure- 
ments of the target of interest, 
using thermometers, thermistors, or 
thermocouples. However, use of 
these devices may, in fact, alter the 
thermal nature of the surface, re- 
sulting in questionable temperature 
measurements. Also, the number 
and frequency of measurements can 
h* prohibitive for large-area phe- 
.omena. An alternative to direct 
measurement of target temperature 
is the use of a radiometer designed 
to measure electromagnetic energy 
in the thermal 1R region. Based on 
certain known or determinable pro- 
perties of the target, measured 
energy can be converted to tempera- 
ture through the use of established 
equations. 

In FY83, the following activities 
were accomplished: 

• A report summarizing the 
underlying theory, as well as 
providing a detailed bibliogra- 
phy dealing with mid/thermal 
IR research, was completed 
(Anderson 1983). The FY83/ 
FY84 experiment design is 
contained in this report. 

so 


• Lowtran-S, an atmospheric 
model developed by the Air 
Force Geophysics Laboratory, 
was implemented on the 
NSTL/ERL mainframe com- 
puter, thus providing the proc- 
essing capability for estimating 
atmospheric parameters re- 
quired for this investigation. 

• Thermal channel ratioing soft- 
ware (post-atmospheric correc- 
tions) was generated, as well 
as software which permits the 
calculation of effective radiant 
temperature from the Thermal 
Infrared Multispectral Scanner 
airborne data collection sys- 
tem. 

• A study site was selected in 
the vicinity of NSTL. Acquisi- 
tion of g r ound truth data, us- 
ing PRT-5 radiometers, was 
scheduled in conjunction with 
a TIMS overflight in the late 
fourth quarter FY83. 

FY84 activities will at least con- 
tinue the thermal portion of the 
study: 

• Initial analysis of 1983 TIMS 
data from the study siu will 
be completed. The results jf 
the analysis will demonstrate 
the applicability of the band- 
ratioing model to a homo- 
geneous land cover. 

• The model will be verified 
using TIMS data to be col- 
lected in 1984 covering a less 
homogeneous land cover type 
(e.g., pine plantation). 


VEGETATION STRESS 
ASSESSMENT 

This RTOP investigation was in- 
itiated to detect plant stress in a 
situation in which some of the vari- 


ables relating to stress in a pine 
plantation can be identified and 
measured. 

The additional data available in 
the mid-lR region of the spectrum 
(I.O-I.3 and 1. 55-2.35 micrometers) 
from TM and other experimental 
scanners require that research be 
conducted to determine the contri- 
bution these regions of the spectrum 
will have on the analysis of vegetat- 
ed areas. In some non-natural vege- 
tated areas such as pine plantations, 
many biological variables (uncer- 
tainties) are removed by the uni- 
formity of planting and by mainte- 
nance procedures. 

The detection of plant stress is of 
importance since, in some cases, it 
is possible to remedy the stressed 
condition. This RTOP examines 
areas within an experimental plant- 
ing of pine near NSTL where there 
are positive indications of differing 
degrees of stress. Within these 
areas, planted in 1960, differing 
amounts of fertilizer were applied 
the first season after outplanting, 
and some plots were not fertilized 
at all. Although 22 years have 
elapsed since this single fertilization 
treatment, certain fertilized plots ex- 
hibit more than twice the wood 
volume of the non-fertilized areas. 
Several strong indications link this 
wood volume increase to differences 
in mycorrhizai involvement between 
the fertilized and non-fertilized 
plots. Ample evidence suggests that 
this type of stress can be detected 
through remote sensing. 

Thematic Mapper Simulator data 
were acquired in December 1982 
and March and May 1933. Primary 
emphasis, however, has centered on 
the acquisition of field information, 
since this must be acquired over a 
relatively long time base to be sta- 
tistically valid. 
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Wood volume data for each tree 
in each plot were calculated, as 
were the average wood volumes per 
plot (each plot is 120 feet square). 

In all areas, * wood volume of 
the once-fertiiized plot is about 
twice that of the non-fertilized plot. 
There is strong evidence which sug- 
gests that this continuing increased 
growth observed in the fertilized 
plot is the result of an increased 
mycorrhizal involvement. Field 
studies to date substantiate this con- 
cept. Mycorrhiza is a symbiotic 
relationship between higher plants 
and various fungi. It is mutualistic 
in that neither organism (pine or 
fungi) can exist without the other. 
There may, however, be differences 
in the degree, extent, and efficiency 
of this partnership. In cases where 
the relationship is not fully exploit- 
ed, the higher plant partner will be 


under more stress than if this were 
not the situation. 

The fungi which are involved in 
this partnership may not be evident 
for much of the ti>'.<e. Under 
favorable conditions, the fungi 
reproduce by producing sporo- 
phores which can be seen in the 
vicinity of the tree with which this 
fungus is a partner. Weekly inven- 
tories were taken during favorable 
weather in the plots, and the 
10-foot by 10-foot square formed 
by each four of the row-planted 
trees was used as the study’s inven- 
tory unit. 

Two such inventory sheets, 
summed up over the entire season 
of appearance, are shown in Figure 
5. The groupings of numbers within 
adjacent blocks suggest the extent 
of the mycelium of the fungus (this 
is the usually underground vegeta- 


tive phase of the fungal organism). 
Note that with both Laccarfa Lac- 
cata and Cortinarlus semisanguin- 
eus , there appears to be a tendency 
for each mycelium to be exclusive 
of the other. This suggests that, 
even within given plots, it may be 
possible to relate individual sporo- 
phore counts for each 10 x 10 
quadrant to determine what correla- 
tions may exist between these field- 
derived data. These correlations will 
be checked in FY84 against the 
TMS data and will be used to deter- 
mine which trees should be sampled 
with a portable field spectrometer. 

Table 6 suggests differences in 
both mycorrhizal spedation and 
extent between the fertilized (high) 
and non-fertilized (low) plots. For 
example, Cortinarius semisanguin - 
eus is much more prevalent in the 
high plots (633 specimens) than in 


L »cc aria laccats 




W2 

W3 

W4 

wT 

Wt 

wr 

wg 

W9 

W10 


HM2 

WTj | 

Most* 









m 



□ 


M2 

2 


1 






a 



B 


Hi 












■ 


H4 

2 











B 


Hi 

2 

1 

■ 

■ 



■ 

B 




z 


N* 

12 

4 

1 



□ 

D 





B 


Hr 



2 











Ml 

f 






D 

B 


2 




Ml 

□ 




D 

■ 

H 

_ 




B 


NiO 

4 





■ 

D 

B 

B 


1 

B 

B 

MU 

1 






B 

B 


t 



B 

HU 

1 



1 

4 


■ 

B 


n 




H 1 1 








L 

C 


□ 





TOTAL: 160 


Cortinarius samisangulnaua 


ksaa 

*2 

W3 

W« 

El 


W7 

Wl 

wt 

MHO 

Wii 

MM? 

W13 

Ml 



■ 

B 


B 

B 


z 

Z 

B 

n 


M2 



j 



1 



B 

B 

B 

Z 


M3 






s 



14 



B 

B 

H4 


B 

B 



D 

B 





_ 


Hi 


B 

B 

B 

4 

B 

B 

B 

B 

B 

B 

B 


Nf> 



B 

B 

B 

s 

9 

B 

B 

■ 

B 

□ 


V 




2 



3 

B 

B 

B 

B 

B 


Ml 






■ 

B 

B 

B 

B 


_ 


Ml 


B 

B 

B 


i 

B 

B 

B 

B 

B 

B 


MiO 

7 

B 

D 

Z 



• 

r 

B 

B 

B 

B 




B 

■ 

B 



B 

B 

B 

B 

B 

Z 



B 

B 





B 

B 

B 

B 

B 

B 



L 

— 


— 



Z 

Z 

_ 

□ 

z 

C 



TOTAL: 653 


Figure 5. Sample Sporophore Inventory Sheets 


51 








APPLIED RESEARCH AND DATA ANALYSIS 


the low plots (104 specimens), while 
Suillus decipitns is much more 
abundant in the low (1,103) than in 


the high (148). These data for the 
first time suggest that C. semisan- 
guintus may be more efficient for 


the tree in the mycorrhizal symbi- 
osis than S. decipiens. 


TABLE 6 

Fertilization Plot Study Area, Harrison Experimental Forest 
(Loblolly Pine Plots) 



BLOCK IV 
PLOT 1 (LOW) 

BLOCK IV 
PLOT 4 (HIGH) 

SPECIES 

Summer 

Winter 

Summer 

Winter 

Amanita flavoconia 

62 

— 

31 



A. rubascans 

307 

— 

66 


Boletus adults 

5 

— 

1 

— 

Qyroporus castanaus 

8 

— 

— 

— 

Tylopilus (Indeclsus) 

1 

— 

— 


T. ptumbaoviolaceus 

2 

— 

— 

— 

Cantharellus cibarius 

26 

— 

— 

— 

C. Clnnabarinus 

24 

— 

— 

— 

Cortinarius cfnnamomaus 

— 

15 

16 

36 

C. samisanguinaus 

— 

104 

— 

653 

Laccaria laccata 

— 

323 

— 

160 

Ttalephora terrastris 

5 

— 

— 

— 

Russula ( amatica ) 

— 

1 

— 

1 

Suillus cothumatus 

— 

2 

— 

2 

Suillus dacipians 

1,103 

13 

148 

— 

Amanita vaginata 

- — 

— 

2 

— 

Amanita vaginata var. alba 

— 

— 

3 

— 

Bolatus projactalius 

— 

— 

6 

— 

Ccllybia maculsts 

— 

— 

— 

6 

Hydnum imbricatum 

— 

— 

— 

1 

Rhizopogon sp. 

— 

— 

— 

2 

Trtcholoma flavovirans 


— 

— 
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JOINT RESEARCH PROJECTS 


Joint research projects represent 
a NASA program area within which 
applied research is conducted jointly 
with selected non-NASA user com- 
munities who have a requirement to 
incorporate remote sensing-derived 
information into their decision- 
making processes. Research projects 
are selected through announcements 
in the Commerce Business Daily 
and subsequent proposals submitted 
by interested organizations. The 
major criteria for selection include: 
(1) the research requirements must 
be of mutual interest to NASA as 
well as the non-NASA organization, 
and (2) the organization cooperating 
in the research must be representa- 
tive of a broader community of 
users. 

During FY82, four joint research 
projects were initiated at NSTL/ 
ERL; the following pages present 
the FY83 accomplishments for each 
project. 


CROP MENSURATION 
AND MAPPING 

The Crop Mensuration and Map- 
ping joint research project is a co- 
operative effort between NSTL/ 
ERL and International Harvester as 
a representative of the farm imple- 
ment industry. The objective of the 
project is to evaluate Landsat The- 
matic Mapper data for mapping 
crop types and estimating crop 
acreages. 

The first year of this project 


(FY82) addressed the development 
of a baseline surface cover inven- 
tory with Landsat MSS data to 
compare with TM-derived inventory 
results. Work completed during 
FY83 focused on the capabilities of 
the Thematic Mapper to discrim- 
inate between crop types in Poinsett 
County, Arkansas, and the feasibili- 
ty of reducing the quantity of data 
analyzed while maintaining high 
mapping accuracies. 

Two methods of spectral signa- 
ture development were employed 
during initial TM analysis: un- 
supervised and supervised. The 
ELAS module SEARCH was used 
to generate unsupervised spectral 
statistics over the Powers Slough 
and Otwell quadrangles study site 
by sliding a 3- x 3-pixel window 
through the data. With the higher 
resolution of TM data, statistics 
derived in this manner are based on 
a sampling unit of 0.8 hectare (2.0 
acres) as compared to 4.0 hectares 
(10.0 acres) when the same tech- 
nique is applied to MSS data. Fifty- 
seven spectral signatures were devel- 
oped over the test site, which cor- 
related with the following seven 
land cover categories: 

• Rice 

• Soybeans 

• Fallow 

• Hardwood 

• Emergent Vegetation 

• Water 

• Clouds 

The resulting classification was 
assessed for accuracy with a file of 


verified ground truth polygons. 
Target crops of rice and soybeans 
were classified correctly at levels of 
98.50 and 98.08 percent, respective- 
ly. 

Supervised spectral signatures 
were developed with ELAS module 
SUPE and the class-naming, ground 
truth data file. Spectral statistics 
developed with this technique were 
used to classify the entire data set 
and resulted in an overall mapping 
accuracy of 92.49 percent correct. 
The target crops, rice and soybeans, 
were classified at 96.64 and 96.30 
percent accuracy levels, respectively. 

Additionally, the supervised sta- 
tistics developed from the seven 
channels of TM data were used to 
examine the utility of each subset of 
channels. The purpose of subset 
selection is to determine the best 
channel combination for any given 
application of the data and reduce 
the number of channels required for 
analysis. The most desirable subset 
of channels was determined to be 
that subset exhibiting the least 
reduction in average transformed 
divergence (ATD) from the original 
channels analyzed. Initially, the 
subset selection routine was applied 
to all seven channels. During this 
analysis, Channel 6 (thermal 1R) 
was found to heavily weight the 
ATD of each subset. This weighting 
was determined to be a function of 
random noise in the data and was 
not due to spectral heterogeneity. 
Thus, Channel 6 was deleted for the 
remainder of the analysis. The over- 
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all ATD for six channels was com- 
puted at 0.25091. The best five- 
channel subset was obtained by 
deleting Channel 3, causing only a 
5.72 percent reduction in ATD. A 
reduction in ATD of 11.19 percent 
was shown by deleting Channels 
and 3 to yield the best overall four- 
channel subset— 2, 4, 5, and 7. 

The assessment of each subset of 
TM channels for discriminating 
crops will be completed during 
FY84 by conducting both supervised 
and unsupervised classifications 
with each subset. The resulting 
classification accuracy figures for 
each subset will then be compara- 
tively analyzed with a test for vari- 
ance to determine the significance 
of each subset test result. 


TIMBER RESOURCE INVENTORY 
AND MONITORING 

The objective of the Timber Re- 
source Inventory and Monitoring 
joint research project with Interna- 
tional Paper Company (IP) is to 
determine the information content 
of Landsat Thematic Mapper data 
for mapping forest stands and iden- 
tifying silvicultural activities. A 
forest stand is defined as “an in- 
dividual body of timber occupying a 
specific area and sufficiently uni- 
form in species composition, age, 
site, and degree of stocking.” In 
addition, a stand must be “separa- 
ble from other areas in the forest.” 
Silvicultural activities include plant- 
ing and restocking, precommercial 
thinning, harvest, site preparation, 
and controlled burning. 

International Paper Company is 
a large-scale forest products indus- 
try operating in North America, 


where it manages some 8.4 million 
acres of forest land. Their current 
inventory is updated on a three-year 
cycle using aerial photography and 
field surveys, which result in detail- 
ed forest data bases. For the pur- 
poses of this project, a study area 
in Baldwin County, Alabama, was 
selected as a representative area of 
southern forest types to take advan- 
tage of IP’s in situ field data. The 
IP data base for the study area in- 
cludes 240 descriptive parameters on 
forest stands as well as silvicultural 
activities. 

Past success with Landsat MSS 
data for industrial forest applica- 
tions has been limited to a general 
forest association level. Even the 
separation of deciduous from conif- 
erous forest has required the anal- 
ysis of seasonal data sets. Overall, 
Landsat MSS data have been con- 
sidered too coarse to provide the 
type of inventory information re- 
quired by industrial forest manag- 
ers. 

During the first phase of this 
joint proiect, an MSS classification 
was developed using seasonal data 
sets to establish a baseline against 
which the TM classification results 
will be compared. While the overall 
accuracy of the MSS classification 
was verified to be 79.2% correct, 
the level of forest information de- 
rived from MSS data was limited to 
general associations, such as mixed 
coniferous and deciduous forest 
communities, and basically substan- 
tiated earlier findings. 

The purpose of the second phase 
of the project is to test existing pat- 
tern recognition techniques using 
Thematic Mapper Simulator data. 
During FY83, several techniques 
were applied to a winter TMS data 


set acquired in February 1982. 
Supervised, unsupervised, and band 
ratioing analysis procedures were 
tested, with the supervised signature 
development technique providing 
the best overall results. The classifi- 
cation resulted in the separation of 
even-aged, single-composition stands 
of loblolly, slash, and long leaf 
pines. Atlantic white cedar, occupy- 
ing the very low areas along major 
drainages, was also discriminated 
from other bottomland hardwood 
stands. 

The verification process employ- 
ed for the accuracy assessment used 
IP timber stand maps, color infra- 
red photography, and field informa- 
tion to identify additional independ- 
ent polygons representative of the 
forest information of interest. These 
polygons were input as a separate 
data file and automatically com- 
pared to the classification data file. 
The overall accuracy for the super- 
vised classification is 88.5%. A 
matrix displaying the interclass 
classification success is presented in 
Table 7. 

The results presented are encour- 
aging in that the TMS data are 
capable of delineating the major 
merchantable timber species of the 
study area and thus address an im- 
portant component of the forest 
stand. These results are made even 
more interesting in that only a 
single winter TMS data set was used 
in analysis. During FY84, a Landsat 
TM data set acquired in late Janu- 
ary 1983 will be used to verify the 
TMS results and to develop new in- 
formation extraction techniques or 
combinations of analysis proce- 
dures, if necessary, based on IP’s 
evaluation of the results thus far. 
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WETLANDS PRODUCTIVE 
CAPACITY MODELING 

The four-year Wetlands Produc- 
tive Capacity Modeling joint re- 
search project, initiated in FY82, 
encompasses the research and 
development of Thematic Mapper 
satellite technology to determine the 
value of wetlands to living marine 
resources. The project is a coopera- 
tive endeavor between the National 
Oceanic and Atmospheric Adminis- 
tration (NOAA) National Marine 
Fisheries Service (NMFS), through 
NMFS’s Southeast Fisheries Center 
(SEFC), and NASA/NSTL/FRL. 

The remote-sensing-based Wet- 
lands Productive Capacity (PC) 
Model (Figure 6) is being refined, 
tested, and analyzed to determine if 
derived results can provide quantita- 
tive information about the value of 
wetland habitats to the production 
of estuarine-dependent fishery 
resources. TM technology is being 
directed toward primary productivi- 
ty assessment and detrital export 
estimation, as interrelated variables 
influencing the determination of the 
trophic values of wetlands. The 
Calcasieu Lake Basin in Louisiana 
is the selected study area for the 
model development. The resulting 
model will interface with another 
model being developed by SEFC to 
describe the energy flow of an estu- 
ary and relate detritus exported 
from the wetlands to the standing 
stocks of juvenile fish and shrimp. 
This NMFS estuarine model will in- 
volve a shrimp productivity index 
and a detritus availability index. 

Because Landsat 4 TM data have 
not been acquired for the study 
area, emphasis was placed on (1) 
testing of the PC algorithms with 
Landsat MSS data, and (2) analysis 
of Thematic Mapper Simulator data 


that had been acquired by aircraft 
during September 1982 and May 
1983. The PC algorithms exist as 
computer programs included in 
ELAS and permit the following 
measurements to be derived from 
digital remotely sensed data (pro- 
gram acronym is indicated in 
parentheses): 

• Surface cover classification 
(MAXL) 

• Water body identification 
(WBOD) 

• Shoreline identi^cation and 
length (SI IN) 

• Distance-' o * ater (DIST) 

The prelL rnary analysis of the 

September ~ i3 data was presented 
at IGA R3£ 83 : n San Francisco in 
August 1983 in a paper entitled 
“Assessment of Wetland Productiv- 
ity Capacity from a Remote-Sens- 
ing-Based Model— A NASA/NMFS 
Joint Research Project” (Butera 
and Frick, 1983). The most impor- 
tant preliminary findings were as 
follows: 

• Regression analyses of TMS 
data showed a stronger rela- 
tionship with percent vegeta- 
tion cover (living and dead' 
than with biomass for esti- 


mating the source of organic 
detritus. 

• Land/water discrimination, 
fundamental to the measure- 
ment of the hydrographic vari- 
ables, was accomplished more 
accurately with TMS than with 
MSS data as compared to a 
standard established through 
aerial photo interpretation. 

FY84 analysis of TMS data will 
focus on the verification of these 
findings, using the TMS data ac- 
quired in May 1983, and also on 
use of band ratioing and other ana- 
lytical techniques. 

NMFS personnel have continued 
analyses to understand some of the 
relationships in their Estuarine Eco- 
system Model. As reported in the 
previously cited IGARSS ’83 paper, 
initial results indicate that salinity 
patterns and estuarine animal inter- 
actions affect the utilization of 
detritus. 

Future efforts will concentrate on 
the calibration of the Productive 
Capacity Model. It will be calibrat- 
ed by multiple regression analysis, 
using organic carbon export as the 
dependent response variable and 
primary productivity, shoreline den- 


TABLE 7 

Mapping Accuracy for Supervised TMS Classification* 


GROUND- 

TRUTH 

CLASSES 

UNCLAS- 

SIFIED 

0 

LOBLOLLY 

PINE 

1 

SLASH 

PINE 

2 

YOUNG 
REGEN. PINE 
3 

LONG LEAF 
PINE 
4 

ATLANTIC 

CEDAR 

6 

1 

0.00 

94.65 

021 

0.00 

4.95 

0.00 

2 

0.25 

0.88 

90.20 

0.50 

8.17 

0.00 

3 

20.62 

0.42 

0.28 

68.62 

12.06 

0.00 

4 

0.25 

4.71 

1.49 

1.49 

92.06 

0.00 

5 

0.33 

0.22 

0.00 

0.00 

0.11 

99.38 


•19 STATISTICS, THRESHOLD 99.9 PERCENT CORRECT OVERALL-88.57 
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Figure 6. Wetlands Productive Capacity Model 


sity, distance to water, and water 
body type as the independent 
variables. Field collection for the 
response variable has been con- 
ducted in FY83 to obtain seasonal 
concentrations of organic carbon in 
the water column for both incoming 
and outgoing tides. The fractions of 
organic carbon sampled and mea- 
sured included dissolved, particu- 
late, and macroscopic detrital 
forms. Dissolved organic carbon 
(DOC) and particulate organic car- 
bon (POC) were derived by auto- 
matic carbon analysis of water bot- 
tle samples. Detrital fragments were 
sampled with nets and dried and 
weighed for calculation of carbon 
content. 

The carbon measurements will be 
defined as organic carbon export by 
adjusting for discharge from four 
hydrologically discrete locations in 


the watershed. The organic carbon 
was sampled and water level and 
current velocity were continuously 
recorded at these four sites. The 
result of calibration, then, will pro- 
vide predictive equations for organ- 
ic carbon export from a typical 
Gulf Coast marsh for four seasons. 
Once the PC model is calibrated, 
the output data will be used as in- 
put data for developing the Estu- 
arine Ecosystem Model by simula- 
tion methods. 


CROP CONDITION ASSESSMENT 
AND MONITORING 

Monsanto Agricultural Products 
Company and NASA/NSTL/ERL 
are participating in a joint research 
project to investigate the capacity of 
Landsat Thematic Mapper data to 


provide agricultural information of 
utility to the agricultural chemicals 
industry. Specific objectives include 
the examination of TM data for 
detecting variations in weed infesta- 
tion levels during the growing sea- 
son and variations in crop residue 
levels in fields between growing 
seasons. 

Work during FY83 focused on 
the examination of weed infestation 
levels in soybean fields in a Gentry 
County, Missouri, study area. The- 
matic Mapper Simulator data show- 
ed that spectral responses of all 
channels except Channel 4 decrease 
with increasing levels of observed 
weed infestation in soybean Fields 
(Figure 7). Although, in general, 
fields with severe weed infestation 
have higher than average cover val- 
ues (percentage of the field covered 
by crop or weeds), the same trends 
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exist in the data for fields of equal 
cover values over 40% cover (Table 
8). It is hypothesized that differ- 
ences in chlorophyll concentrations 
and plant water content between 
vigorously growing weeds and the 
soybean crop account for the 
changes in detector response as 
weed infestation increases. 

Using simulated TM data, three 
classification methods were examin- 
ed for their effectiveness in deter- 
mining broad levels of weed infesta- 
tion in soybean fields. The first 
method was a supervised classifica- 
tion scheme in which classification 
statistics were derived by merging 
statistics from fields of similar weed 
infestation levels. Through this 
method, three sets of statistics rep- 


resenting severe, moderate-light, 
and no weed infestation were cre- 
ated. The statistics were then used 
to classify an independent set of 
verification fields of known infesta- 
tion levels. The classification was 
59.3% correct in identifying fields 
with the majority (greater than 
50%) of their pixels in the correct 
weed infestation level according to 
ground truth data. 

The second method consisted of 
developing three to four statistics 
per weed infestation level by search- 
ing each level with an unsupervised 
pixel clustering algorithm. The re- 
sulting statistics were used to classi- 
fy the fields. Using this method, 
classification accuracy was increased 
to 66.7%. The increase in accuracy 


was due to the use of classification 
statistics which better approximated 
the actual centers of variation 
found within the weed infestation 
levels. The statistics indicate a het- 
erogeneity in soybean fields with 
light-moderate and severe weed in- 
festation. Even within these fields, 
there are centers of variation which 
represent pixels of pure crop, as 
evidenced by the similarity to the 
statistics from fields with no weed 
infestation. 

The third method consisted of 
assigning the statistics developed in 
the second method to one of three 
weed infestation levels based on 
their statistical similarity to mean 
signatures for each weed infestation 
level. The resulting statistics when 
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Figure 7. Spectral Signatures of Soybeans Within Four 
Weed Infestation Levels 


59 



JOINT RESEARCH PROJECTS 


used to classify the verification 
fields yielded an accuracy of 70.3%. 
The increase in accuracy over the 
first two methods results from a 
decrease in the standard deviations 
of the statistics representing the 


weed infestation levels. 

The Crop Condition Assessment 
and Monitoring joint research proj- 
ect will continue in FY84. Projected 
activities include continuation of the 
weed infestation work, using corn 


as the target crop, and an examina- 
tion of the capability of simulated 
TM data to detect variations in 
crop residue levels with a late Oc- 
tober data set. 


TABLE 8 

Mean Channel Response for Weed Infestation Levels within Four Coverage Levels in Soybean Fields 


COVERAGE 
VALUE (%) 

WEED 

INFESTATION 

LEVEL 

TMS CHANNELS 

1 

2 

3 

4 

5 

6 

7 

80-100 

SEVERE 

102.81 

113.50 

91.57 

147.35 

142.21 

90.77 

106.08 


MODERATE 

103.42 

114.58 

91.77 

162.12 

148.01 

91.06 

103.22 


LIGHT 

105.61 

112.11 

98.17 

147.11 

157.07 

93.95 

117.29 

60-80 

SEVERE 

103.18 

1 1 1 .55 

101.18 

110.43 

144.64 

93.62 

117.34 


MODERATE 

107.39 

117.64 

110.69 

115.13 

155.63 

97.25 

139.49 


LIGHT 

107.02 

116.30 

105.74 

124.66 

155.19 

96.11 

137.60 


NONE 

110.97 

122.84 

123.36 

110.84 

174.22 

101.67 

145.32 

40-60 

MODERATE 

103.59 

111.19 

99.56 

111.68 

142.67 

94.07 

136.32 


LIGHT 

109.37 

120.50 

112.26 

124.66 

158.82 

97.42 

128.34 


NONE 

111.98 

129.84 

132.56 

121.64 

185.56 

102.85 

151.39 

0-40 

MODERATE 

111.37 

123.46 

120.99 

115.86 

163.28 

99.51 

136.28 


LIGHT 

113.87 

128.02 

129.02 

123.07 

171 .56 

101.45 

146.62 


NONE 

108 55 

118.03 

116.98 

103.98 

164.05 

99.62 

174.89 
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Four permanent test areas within 
the NSTL/ERL region have been 
identified for development under 
the Southern Test and Evaluation 
Program. The first area to be devel- 
oped is the Gulf Coastal Plain (Fig- 
ure 8), within which several project 
activities have been initiated. 

The Gulf Coastal Plain is charac- 
terized by a wide variety of environ- 
mental and cultural conditions, 
which led to its selection as a test- 
ing ground for remote sensing tech- 
nology. During Fiscal Years 1982 
and 1983, work was conducted on a 
specific test site in Baldwin and 
Mobile Counties, Alabama, which 
is representative of the Gulf Coastal 
Plain physiographic region. Activi- 
ties completed during this time in- 
cluded Landsat MSS-derived land 
cover inventories using unsupervised 
spectral signature development tech- 


niques for: (1) Mobile/Baldwin 
Counties, (2) Baldwin County, (3) 
Baldwin County Coastal Zone, and 
(4) the City of Mobile. In addition, 
a standardized data base for Bald- 
win County was constructed and, 
together with the MSS-derived land 
cover inventories, will be used to 
test and evaluate emerging sensor 
technology and new data analysis 
algorithms, including data base 
models. 


CORRIDOR ANALYSIS 

During FY83, the Baldwin Coun- 
ty data base was used to test a 
weighted simulation model develop- 
ed at NSTL/ERL from a software 
package developed by Environmen- 
tal Systems Research Institute 
(ESRI). The model was composed 


of five data planes geographically 
referenced to the UTM coordinate 
system on a 30-meter grid cell size 
and used to determine relative levels 
of trafficability for each cell. Data 
planes used to model the area in- 
cluded: 

• Landsat MSS-derived land 
cover 

• U.S. Soil Conservation Service 
(SCS) detailed soils survey 

• National Cartographic Infor- 
mation Center (NCIC) digital 
topographic data 

• Digitized primary and secon- 
dary roadways 

• Digitized urban centers 

The original data comprising each 
data plane were very detailed and 
contained more information than 
was required for the application. In 
order to incorporate these data into 
the weighted simulation model, the 
variables within each data plane 
were aggregated into relevant cate- 
gories based on hypothetical criteria 
for determining trafficability. 

In some cases, this process entail- 
ed aggregating the detailed informa- 
tion contained in the original data 
plane. For example, 100 different 
soil types were combined into 6 
load-bearing strength categories, 
and the 38 spectral classes derived 
from Landsat MSS data were com- 
bined into 9 land cover classes. In 
other cases, additional information 
was derived or computed from the 
original, as in the elevation data. 
From the NCIC digital terrain data, 
elevation is the original data plane 
variable from which aspect and 
slope can be derived. For the appli- 
cation, percent of slope was com- 
puted as a limiting factor for traf- 
ficability. The importance of trans- 



Figure 8. Gulf Coastal Plain Test Area 
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TEST AND EVALUATION 

portation networks and urban cen- 
ters was a function of distance, 
which was computed using an algo- 
rithm which determines the Euclid- 
ian distance between specified units. 

In weighted simulation modeling, 
each ceil representing the categories 
within each data pidne was assigned 
a ranking value (X value) between a 
range from 0 to 10, with 0 repre- 
senting the best condition and 10 
representing the worst, as depic**d 
in Figure 9. The conditions placed 
on corridor development in this ex- 
ercise were assigned arbitrarily. For 
example, within the land cover data 
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plane, 38 spectral classes were com- 
bined to form 9 categories. These 
were ranked in terms of trafficabili- 
ty, with those categories offering 
the least resistance to corridor de- 
velopment (cultivated lands, pas- 
ture, etc.) assigned the lower values. 
In each case within the remaining 
data planes, the lower ranking 
values were assigned to the areas 
offering the least resistance to cor- 
ridor development. 

A significance value (N-weight) 
was then assigned to each data 
plane. The N-weight was multiplied 
by the X value for each cell (ex- 


pressed as Xj Nj * acceptability 
value) to determine the relative traf- 
flcability value within each data 
plane. The data planes were over- 
laid and the corresponding cell ac- 
ceptability values within each data 
plane were then summed to deter- 
mine a range of trafficability 
(Figure 9). 

For the actual corridor selection, 
given Point A as the departure 
point and Point B as the target, the 
algorithm defines a minimum dis- 
tance (path of least resistance) from 
Point A to Point B (Plate 9). In 
Plate 9, all data planes had differ- 


(1) CATEGORIES AND COR- 
RESPONDING 50-METER 
CELLS ARE ASSIGNED AN X 
VALUE 


EXAMPLE 


(2) X VALUE THEN MULTIPLIED 
BY WEIGHT (N VALUE) 

ASSIGNED TO CORRE- 
SPONDING DATA PLANE 
(X,N,) + - <X n N n ) * ACCEPTABILITY VALUE 
WEIGHTED DATA PLANES 


(3) RESULTS SUMMED TO 
DEFINE RANGE OF 
ACCEPTABILITY 


RESULTS 
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Figure 9. Weighted Simulation Model 
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ent weights of importance assigned. 
Note that only when the urban ex* 
elusion zone (light blue) was assign- 
ed a 30^c weight did the corridor 
selection algorithm 'void the urban 
center. A 10V# weigt.i did not result 
in avoidance of the urban area. 

The purpose of this hypothetical 
example was to test weighted simu- 
lation modeling within the ELAS 
software package. The significance 
lies in the flexibility of ELAS: 

(1) to handle data from different 
sources available at different scales, 

(2) to standardize by geographically 
referencing each data set to the 
UTM coordinate system, and (3) to 
overlay each data plane for analy- 
sis. No restriction is placed on the 
number or type of data planes 
which can be used. 

As decision-makers place increas- 
ing reliance on automated ap- 
proaches for synthesizing more and 
more data, weighted simulation 
modeling represents a valuable tool 
for developing a range of options 
for consideration. Within the Gulf 
Coastal Plain test area, this model 
will be applied during FY84 to 
assessing coastal land capability by 
identifying areas for preservation, 
conservation, or utilization based 
on a number of variables ranging 
from environmental to political 
considerations. 


URBAN STUDIES 

Since the launch of the first 
Landsat satellite, geographers and 
land use planners have been inter- 
ested in using space-acquired data 
for urban mapping. A diverse range 
of data analysis techniques has been 
applied to Landsat MSS data by 
researchers, ranging from tradition- 
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al photo interpretation to digital im- 
age processing algorithms, with 
varying degrees of accuracy. The 
limitations of MSS data content for 
delineation of urban morphology 
have been due to the rather coarse 
spatial and spectral resolutions of 
the data. 

The advent of the Thematic 
Mapper, with its improved spectral, 
spatial, and radiometric resolution, 
holds promise for those concerned 
with mapping urban milieu, par- 
ticularly with respect to understand- 
ing the impact of urban growth as 
it changes the character of the sur- 
rounding environment. To support 
the evaluation of TM data for 
urban area analysis, several studies 
were initiated during FY82 and 
FY83. Within the Gulf Coastal 
Plain test area in Alabama, Mobile 
is currently providing a permanent 
site to test advanced sensors. Com- 
plementary to the Gulf Coast work, 
NSTL/ERL has been cooperating 
with the University of Georgia and 
the University of New Orleans in an 
effort to quantify the capability of 
TM data for urban studies. 

Under a research grant, the Uni- 
versity of Georgia’s Department of 
Geography has examined the utility 
of Thematic Mapper Simulator data 
for discriminating urban features in 
the Atlanta, Georgia, metropolitan 
area. TMS image data have been 
evaluated to assess the potential of 
30-meter resolution multispectral 
data for mapping land cover. Visual 
photo interpretation of the data 
permitted 33 land covers to be map- 
ped at Level 11 and Level III (urban 
classes only) from the Anderson, et 
al., USGS classification system. A 
point sample verification of the 
da..j, predicated on comparison 
with aerial photography, has indi- 
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cated that the TMS-derived land 
cover map has an overall accuracy 
of 80^» correct classification. 

Digital classifications were also 
produced from the TMS data 
using supervised techniques. An 
assessment of the optimal channels 
for classification of urban areas was 
determined by employing average 
transformed divergence and prin- 
cipal components analysis. Channels 
2, 4, and 7 were determined to be 
the best three-channel data set for 
supervised classifications; the prin- 
cipal components approach produc- 
ed a three-component set (principal 
components 1, 2, and 3) which gave 
accuracies comparable to those of 
the Channels 2, 4, and 7 ruw data 
set. 

Additionally, one of the first 
seven-channel TM data sets avail- 
able to investigators from Landsat 4 
was acquired over portions of Ar- 
kansas and Tennessee on August 
22, 1982. In support of NASA's 
LI DQ A program, researchers at 
NSTL/'£RL conducted a prelimi- 
nary evaluation of the TM data for 
discriminating and delineating utban 
features occurring within Union 
City, a small city in northwestern 
Tennessee which has a population 
of approximately 15,000 and en- 
compass about eight square miles 
(20.7 km 3 ) within its city limits. 

A qualitative examination of the 
enhanced photo-interpretive quali- 
ties of TM data was conducted. In 
a false color composite image pro- 
duced for the Union City area using 
TM Channels 2, 3, and 4, the 30- 
meter ground resolution of the data 
readily permits photo interpretation 
of individual structures and other 
urban features. This is a significant 
advancement over what has been 
available with MSS data, where the 
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80-meter spatial resolution allowed 
photo interpretation of only large 
urban phenomena. 

The TM data were also evaluated 
for use in a digital classification of 
Union City. A pixel-to-pi;:el signa- 
ture development algorithm was 
employed to derive statistics for a 
polygon of data delineated around 
the Union City urban area. Thirty- 
nine spectral signatures were derived 
using TM Channels 2 through 7 as 
input data. 

The-e spectral responses were 
then classified and the classes were 
grouped, utilizing ground truth and 
ancillary data, into six specific land 
cover categories: (1) roads and inert 
materials, (2) commercial and in- 
dustrial development, (3) residential 
development, (4) agricultural and 
bare soil, (5) transitional or grass- 
land areas, and (6) forested areas. 


Plate 10 shows the classification of 
Union City as delineated by the 
polygon that was used to define the 
urban area for signature develop- 
ment purposes. 

An accuracy analysis was com- 
puted for the classification based on 
the evaluation ol ground verifica- 
tion polygons with known land 
covers selected from the data. Be- 
cause of the difficulty in defining 
precise ground truth polygons for 
the transitional/grassland and 
roads/inert materials classes, poly- 
gons for these categories were com- 
bined with the agricultural/bare 
soil, ccmmercial/industrial, or resi- 
dential categories. Table 9 gives the 
accuracies for the resulting four ma- 
jor land . pve r types used to esti- 
mate c!as ification accuracy. 

A more in-depth discussion on 
the classification techniques used 


and the results from analysis for 
Union City is presented in Quat- 
trochi, et at., 1982, and Quattrochi, 
1983a and 1983c. This preliminary 
evaluation was presented at the 
Early Results Workshop at God- 
dard Space Flight Center in Febru- 
ary 1983, and at the URISA Con- 
ference held in Atlanta, Georgia, 
during August 1983. 

The Union City investigation was 
an initial attempt at defining the 
utility of TM data for delineating 
discrete urban features. Present 
research at NSTL/ERL focuses on 
the analyst of TM data for dis- 
criminating urban morpho'iog’’ w ! *‘ 
in a large metropolitan are- -obi 1 
has been established as ont 
nent study area for urban ». es..ga- 
tions. 

The TM data used for this re- 
search were collected on October 
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Computed Classification Accuracies— Union City Study Area* 
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PLATE 10 



Computer-implemented land cover classification of Union City. Tennessee, derived from 
Thematic Mapper dat» acquired on August 22. 1982 


JlSCRlMINATION OF URBAN FEATURES WITHIN A SMALL CITY 


i 

■ 



TEST AND EVALUATION 


2 f, 1982, and provide the basis of a 
joint investigation between NSTL/ 
ERL and the Department of Geog- 
raphy at the University of New 
Orleans. To evaluate the digital 
classification products that could be 
derived from the data, a computer- 
implemented land cover classifica- 
tion was produced from TM Chan- 
nels 1, 2, 3, 4, 3, and 7 for the 
Mobile area. Derivation of the clas- 
sification proceeded in t wo phases. 
Spectral signatures were first devel- 
oped by passing a 3- by 3-pixel 
moving window through the data. 
The six channels of raw data were 
then classified using the signatures 
produced from this technique. An 
analysis of the classification showed 
that 18.5?o of the data were left 
unclassified. This unclassified class 
fell predominantly within the built- 
up areas of Mobile. 


The second phase of data analy- 
sis employed a pixel-by-pixel signa- 
ture development program to derive 
spectral statistics for only those pix- 
els which fell within the unclassified 
area of the previous classification. 
These signatures were used to com- 
plete classification of the data. The 
result of this two-phased process 
was a composite classification of 
the Mobile area with 12 land cover 
classes (Plate 11). Although an 
analysis of the classification for 
Mobile is still in progress, an inter- 
im evaluation of the classification 
has been documented in Quattrochi, 
1983b and 1983c; these initial re- 
sults were presented at the ERIM 
Symposium at Ann Arbor, Michi- 
gan, in May 1983 and at the 
URISA Conference in Atlanta, 
Georgia, in August 1983. 

Findings to date indicate there is 


more informational content in the 
TM data than can be extracted 
through existing spectral digital 
analysis techniques, and that the 
spatial attributes of the data have 
yet to be fully addressed. Further- 
more, due to the combined increas- 
ed spectral and spatial resolutions, 
there is more information in the 
resulting classification over Mobile, 
for instance, than can be adequately 
assessed by standard methods for 
computing classification accuracy. 

Work at NSTL/ERL in coopera- 
tion with the University of Georgia 
and the University of New Orleans 
will continue during FY84 to evalu- 
ate TM data for urban studies and 
to develop new information extrac- 
tion techniques which can capture 
both the spectral and spatial infor- 
mation content of the data. 


PRECEDING PAGE BLANK NOT FILMED 





71 


TEST AND EVALUATION 
Urban Studies 


ORIGINAL PAGE 
COLOR PHOTOGRAPH 


PLATE 11 




: - v % 


' > 




— 


R E S3 DEN T I-R L < m£ R T, ; V 
RESlbfeNf*m>F6PE5T 
. UP Eft. I i»fiSSLHrip . : 

£ INDUS T<COHHE 
■E:fiRE" SO 
INERt-O flgfkv 
• FOREST ' • * - • -* / 

. BRUSH LfiTjD- TRRNSI TIjONAL 
CUL TURt ■ORHi^S LAND 
' fRSH LrRHD' ; wVt LflNb 
•WATER 3-3 . 

'1 OTHER LflNCT £'0to£RS* 






• _ 




t / 




This classification of the Mobile metropolitan area, from TM 
data acquired October 27, 1982, was produced in two 
phases— first by using a 3- by 3-pixel moving window to 
develop signatures for “homogeneous" land covers and 
then by employing a pixel-by-pixel signature development 
dlgoriihm which derived signatures for tne residual, spec- 
trally heterogeneous urban area 


TM LAND COVER CLASSIFICATION OF MOBILE, ALABAMA, AND ENVIRONS 
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NSTL has actively participated in 
the preliminary Space Station defi- 
nition activities through working 
groups, short-term technology 
studies, and Program Operating 
Plan (POP) preparation, submis- 
sion, and advocacy. 

In FY83, significant contribu- 
tions were made to the Space Sta- 
tion activities through the following 
working groups: Operations Work- 
ing Group, Mission Requirements 
Working Group, Autonomy Work- 
ing Group, Data Management 
Working Group, and Software 
Working Group. 

NSTL/ERL was asked to partici- 
pate in one of 14 short studies ini- 
tiated by the Space Station Task 
Force Concept Development Group. 
The User Accommodation data sys- 
tem study was completed in four 
weeks by drawing from existing per- 
sonnel experience. The primary pur- 
pose of the study, which was done 
in parallel with other NASA cen- 
ters, was to determine if there were 
any significant cost drivers in the 
user interface portion of the Space 
Station Data Management System. 

A comparison of the individual 
centers’ findings indicated no signif- 
icant cost drivers. 

The NSTL prepared and present- 
ed four Program Operating Plan 
proposals to the Space Station Task 


Force Peer Review Group. The four 
studies were Data Management Sys- 
tem Test Bed; Evolutionary Soft- 
ware Study; Expert System: Change 
Detection System; and Closed Eco- 
logical Life Support System Study. 
Although at this writing no final 
decision has been made on approval 
of any of the proposals, preliminary 
reports are encouraging and NSTL 
anticipates acceptance of one or 
more of the proposals. 

NSTL will continue to actively 
support Space Station related ac- 
tivities and intends to undertake a 
significant role in Space Station 
development once it becomes an ap- 
proved program. 

NSTL/ERL has been selected to 
participate in the Data Management 
System Test Bed with Johnson 
Space Center, Marshall Space Flight 
Center, Kennedy Space Center, Jet 
Propulsion Laboratory, Goddard 
Space Flight Center, Ames Research 
Center, and Langley Research 
Center. NSTL’s effort will be con- 
centrated in the areas of user par- 
ticipation, sensor simulation, data 
and command, file protection, and 
software upgrade techniques. 
NSTL/ERL holds membership on 
the Space Station End-to-End Infor- 
mation System Test Bed Steering 
Committee. 

A summary of NSTL’s Space 


Station involvement dunng FY83 is 
as follows: 

• Committee Participation 

• Operations Working Group 

• Autonomy Working Group 

• Higher Order Language 
Working Group 

• Data Management Working 
Group 

• Mission Requirements 
Working Group 

• Software Working Group 

• Data Management System 
Test Bed Steering Commit- 
tee 

• Concept Development Group 

Short Study 12: User Accom- 
modation Data Systems 

• Assigned to study team 

• Extracted data from person- 
nel experience base 

• Documented study 

• Presented findings to com- 
bined group at NASA 
Headquarters 

• Space Station Task Force 

POP Call 

• Data Management System 
Test Bed 

• Evolutionary Software 
Study Submittal 

• Expert System: Change 
Detection System 

• Closed Ecological Life Sup- 
port System Study 
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W. G. Cibula, 688-1913 
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SourcM of Additional Information 

For additional copies of this 
report, contact: 

E. F. Zetka 

NASA/NSTL Earth Resources 

Laboratory 

NSTL, MS 39529 

(601) 688-3830 

For additional information on 
the topics discussed, contact the 
principal authors listed below. 

SENSORS AND DATA SYSTEMS 

Thermal Infrared Multispectral 
Scanner - G. F. Flanagan, 
688-3588 

Advanced Sensor Design and 
Operation - G. F. Flanagan, 
688-3588 

Software Development and Man- 
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and Monitoring - G. S. 
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plications - K. J. Langran, 
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APPLIED RESEARCH AND 
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geneous Mine Features - 
D. A. Quattrochi, 688-1919 
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search - H. B. Musick, 
688-1918 

Geological Mapping - 
D. L. Rickman, 688-1920 
Scene-to-Scene Registration - 
J. E. Anderson, 688-1909 
Topographic Corrections - 
D. L. Rickman, 688-1920 
TM Preprocessing Technique - 
J. E. Anderson, 688-1909 
Geo botanical Application - 
W. G. Cibula, 688-1913 
Scene-to-Map Registration - 
D. D. Dow, 688-1914 
Microwave Sensor Analysis - 

S. T. Wu, 688-1922 


Crop Mensuration and Map- 
ping - D. P. Brannon, 
688-2043 

Timber Resource Inventory and 
Monitoring - C. L. Hill, 
688-2047 

Wetlands Productive Capacity 
Modeling - D. D. Dow, 
688-1914 

Crop Condition Assessment and 
Monitoring - G. J. Irish, 
688-1907 
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